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A theory  is  presented  for  determining  critical  loads  for  cylindrical 
shells  subjected  to  Large  underwater  ex|)losions.  The  theory  treats 
dvn.imic  elastic-plastic  buckling  with  fluid-structure  interaction  caused 
by  a transverse  incident  pulse  in  water.  The  approach  is  by  modal  analysis 
with  the  displacement  and  velocity  distributions  at  tlie  end  of  the  initial 
elastic  response  forming  the  initial  conditions  of  the  subsetjuenl  plastic  — - 
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20.  ABSTRACT  (Continued) 

‘respon.se.  The  transition  from  elastic  to  plastic  response  is  determined 
by  the  hoop  mode. 

Experimental  and  theoretical  final  deformed  shapes  and  pressure 
histories  are  similar. 

The  response  of  a stiffened  shell  to  a rectangular  pressure  pulse 
is  characterized  by  isodamage  curves  in  the  pressure-impulse  plane. 
Results  show  steep  damage  gradients  when  pressures  and  durations  exceed 
the  values  that  cause  incipient  damage,  which  is  in  agreement  with 
experimental  observation. 
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SUMMARY 


A theory  is  presented  tor  determining  critical  load  curves  for 
cylindrical  shells  subjected  to  large  underwater  explosions.  The 
theory  treats  dynamic  elastic-plastic  buckling  with  fluid-structure 
interaction  caused  by  a transverse  incident  pulse  in  water.  Specific 
results  are  generated  for  a rectangular  incident  pressure  pulse 
(instantaneous  rise,  constant  pressure,  instantaneous  decay)  as  a 
reasonable  idealization  of  a shock  wave  from  a large  underwater 
explosion  with  surface  cutoff. 

Wavelengths  observed  in  experiments  (Fig  3.11)  indicate  that  the 
predominant  ljuckling  mode  forms  during  the  initial  elastic  phase  of 
deformation.  This  mode  enters  the  initial  conditions  of  the  plastic 
phase  of  deformation  and  is  the  preferred  mode  for  further  amplification. 
Tlicrefore,  the  theoretical  approach  is  to  treat  tlie  initial  elastic 
response  of  the  shell  by  model  analysis  until  plasticity  occurs  in  the 
hoop  mode,  and  to  use  the  displacement  and  velocity  distributions  at  this 
time  as  initial  conditions  for  the  subsequent  plastic  response.  Amplifi- 
cation of  the  preferred  buckling  mode  is  assumed  to  continue  during  the 
plastic  phase  until  the  hoop  mode  roaches  the  maximum  inward  displacement. 

Kxpcrimental  and  theoretical  final  deformed  shapes  are  similar 
(as  shown  in  1 iguros  .3.11  and  5.9),  indicating  the  validity  of  the 
a|)proach.  Kxpe ''imen  t al  and  theoretical  pressure  histories  are  also  simi- 
l.ir  (I'igures  3.12,  5.1(5,  A-3). 

Structural  .damage  is  represented  by  an  isodamage  curve  in  the  P-I 
( pressure- impul  se ) plane.  lor  an  incident  rectangular  |uilse  of  prcssuiH' 

P and  duration  f , the  impulse  is  1 P'f . Damage  is  taken  as  the  maximum 


permanent  radial  displacement  ; in  dimensionless  terms  ^ may  be  a trac- 
tion ot  the  shell  radius  or  a multiple  ot  the  shell  thickness.  A pressure 
impulse  curve  is  thereioro  the  curve  in  the  P-I  plane  on  which  is  a 
constant;  the  curve  connects  rectangular  pulses  ol  equal  el lect iveness 
that  range  from  ideal  impulses  to  step  loading.  Pressure-impulse  curves 
(Kigure  6.2)  wore  generated  tor  an  externally  stitlened  cylindi’ical  shell 
(Table  5.1,  Figure  1.2a). 

It  is  shown  by  example  tiiat  the  excess  impulse  is  not  a constant 
along  an  isodamage  curve  but  may  vary  by  as  much  as  a tactor  ol  two  ; 
this  implies  that  e.xcoss  impulse  is  not  as  accurate  a measure  ot  the 
damage  potential  of  underwater  shock  waves  as  the  P-1  representation. 

Plots  ol  the  relationship  between  damage  ^ and  pressure  or  duration 
(Figures  6.3  and  6.1)  show  how  damage  gradients  steepen  with  increasing 
P and  T,  in  conlormance  with  o.xpor imen t al  observation  (Figui’es  3.13  and 
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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement. 


To  Convert  From 

To  i 

Multiply  By  ! 

angstrom 

meters  (m) 

1. 000  000  .X  K -10  i 

atmosphere  (normal)  1 

kilo  pascal  (kPa) 

1 013  25  X E +2  .> 

bar 

kilo  pascal  (kPa) 

1 . 000  000  X E 4 2 

bam  1 

meter^  (m^) 

1. 000  OOO  X E -28 

British  thermal  unit  (thermochemical)  ^ 

1 joule  (J) 

1 1 . 05-1  350  X E 4 3 

calorie  (thermochemical)  1 

1 joule  (J)  i 

j 4. 1S4  000 

1 ' 

; cal  (thermochemical)/cm**  j 

2 2 

mega  joule/m  (MJ/m  ) | 

! 4.  1S4  OOO  X E -2 

I curie  I 

•giga  becquerel  (GBq) 

3.700  OunXF.jl  1 

degree  (angle) 

radian  (rad) 

1 . 745  529  .K  E -2  | 

degree  Fahrenheit 

1 degree  kelvin  (K) 

(t"  f 4 459.07)/!.  8 ‘ 

electron  voU 

joule  (J) 

1.002  19  X E -19  1 

1 1 

! joule  (J) 

1. 000  000  X E -7 

erg/second 

watt  (W) 

1.000  000  X E -7 

foot 

meter  (m) 

3.  04S  000  X E -1 

j foot-pound -force 

1 joule  (J) 

1.355  818 

I gallon  (C.  S.  liquid)  1 

meter^  (m^) 

3.  785  412  X E -3 

j inch  1 

^ meter  (m)  i 

i 2.  540  000  X E -2 

1 jerk  1 

joule  (J)  ' 

1 . 000  000  X E 4 9 

1 jouleAilogram  (J/ltg)  (radiation  dose  i 

] absorbed)  | 

Gray  (Gy) 

1. 000  000 

1 kilotons  j 

terajoules 

4.  183 

! kip  (1000  Ibf)  ! 

newton  (N)  i 

4.  448  222  X E 43 

; kip/inch^  (ksi)  | 

kilo  pascal  (kPa)  ' 

(»  894  757  X E 43 

j ktap  i 

newton -second/m'^ 

(N-s/m^) 

1. 000  000  X E 4 2 

micron  ! 

meter  (m) 

1 000  000  X E -n 

mil  j 

meter  (m)  i 

2.  540  000  X E -5 

mile  (international)  i 

meter  (m)  1 

1 . 009  344  X E 4 3 

{ ounce 

kilogram  (kg)  t 

2. 834  952  X F -2 

5 pound -force  (lbs  avoirdupois)  | 

newton  (N)  ’ 

4.448  222 

pound -force  inch  | 

newton-meter  (N*m) 

1. 129  848  X E -1 

! pound -force/inch  j 

newton/metcr  (N/m)  ' 

1.  751  288  X E 42 

1 pound -force /foot'^  j 

kilo  pascal  (kPa) 

4.  788  020  X E -2 

kilo  pascal  (kPa)  j 

0.  894  757 

i pound-mass  (Ibm  avoirdupois)  j 

kilogram  (kg)  j 

4. 535  924  X E -1 

1 p(>mid-mas.s-foot^  (moment  of  inertia)  ' 

kilogram-meler“  1 

1 : 

(kg-m'*^)  1 

4. 214  Oil  X E -2 

1 3 

pound -mass/Too!  i 

1 kilogram  .''mete r^ 

1 

(kK/m-’i 

1. fiOl  H40  X E 4 1 

rad  (radiatimi  dose  absorbed)  ' 

••Gray  (Gy) 

1. 000  000  X E -2 

j roentgen 

coulomb  Ai  log  ram 

(C/kg) 

2.  579  7(50  X E -4 

shake  l 

second  (s) 

1, 000  000  X E -8 

j 1 

' kih^ram  (kg) 

1. 459  390  X E 4 1 

1 torr  (mm  Hg.  0“  C) 

kilo  pascal  (kPa)  | 

; 1.333  22  X E 

•Tho  bcrquorcl  is  the  SI  unit  of  radioactivity;  1 liq  1 event  ^s. 
••The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 
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1.  INTKomX’TlON 


Uamatio  to  a Bubmeri;cd  submarine  as  the  I'esult  <>l  an  underwatei' 
explosion  is  commonly  divided  into  two  principal  classes,  hull  damat;e 
and  equipment  damage.  Severe  hull  damage  aftects  the  seaworthiness  ol 
the  submarine;  damage  to  equipment,  i.e.,  machinery,  weapons  systems, 
and  other  items  that  are  not  part  ol  the  hull,  afl'ects  mobility  and 
lire  power. 


Damage  to  the  cylindrical  pressure  hull  ol  a submerged  submarine 
by  an  underwater  nticlear  explosion  is  currently  predicted  by  the  hlxcess 
Impulse  Rtile  (described  in  Appendix  A).  This  is  an  empirical  lormula 
for  the  maximum  stiffener  deflection  and  is  based  on  the  correlation  ol 
test  results  lor  a variety  ol  hulls  loaded  by  approximately  rectangtilar 
pressitre  pulses  that  simtilale  those  1 rom  large  c.xplosions  at  large 
standoff.  Test  rcstilts  correlatc'd  by  the  Kxcess  Impulse  Rule  exhibit 
considerable  scatter.  Moref>ver,  becatise  ol  the  empirical  nature  of  the 
Excess  Impulse  Rule,  its  appl  i cabi  1 i t v to  the  prediction  of  tlamage  out- 
side the  range  ol  the  tests  on  which  it  is  based  is  questionable.  In 
particular,  beyond  a certain  depth,  the  Excess  Impulse  Rule  indicates 
that,  lor  a given  pre.ssure  pul.se,  st  i 1 tenner  delleition  de<'  reases  with 
depth  of  submergence.  Since  hv<lrostat  ic  pressure  incrt'ase-s  with  depth, 
it  is  intuit  ivelv  obvioii.s  that  the  depth  dependence  ol  the  Exces.s  Impulse' 
Rule  is  incorrect. 


The  purpose  ot  the  work  reporti'd  hi  n is  pri  vidi  a I heorv  lor 
correlating  available  expi- r tmen  I a 1 results  aiitl  tor  t n 1 1 rpo  1 ,i  i i ng  atid 


extrapolating  the  results  to  other  situations, 
and  extrapolation  valid,  it  is  essenti.il  lh.it 


ions.  lo  ."..he  n t rpol  ,1 1 1 on 
lh.it  the  th,  ,.iv  ’><  li.i'-i'd  on 


r 


1 


I 

the  correct  physical  mechanisms.  To  detei’mine  I lie  mechanisms  on  which 
I ho  theory  should  be  based,  we  examined  hull  ilo 1 orma t i ons  in  scale  models 
produced  bv  underwater  explosions.  From  these  observations  the  theory 
was  developed. 

Critical  load  curves  tor  military  structures  under  attack  loads 
applied  by  t;as  pressure  (air  blast,  lor  example)  have  been  developed 
over  the  past  decade  or  so  in  terms  ol  critical  peak-pressure  impulse 
curves.  The  utility  ol  those  curves  i'oi'  loads  applied  by  p;as  pressure 
led  to  the  sujtKestion  that  similar  critical  load  curves  be  developed 
lor  loads  resulting  irom  underwater  explosions. 

In  structures  loaded  by  ttas  [tressure,  the  compressilti  1 ity  ot  the 
■ias  makes  the  load  essentially  independent  ol  tlie  response.  However, 
loads  resultinti  Irom  sliock  waves  in  water  depend  strons>ly  on  the  response. 

This  leature  — the  f luid-structtire  interact  ion--Kreat  ly  comiilicates  the 
determination  oi  critical  load  curves. 

In  our  first  attempt  to  determine  critical  load  curves  lor  under- 
water explosions,^  we  used  a riiiid-plastic  representation  ol  the  struc- 
ture and  determined  the  critical  loads  at  Ute  lontt-du  ra  t i on  and  short- 
duration  extremes,  taking  approximate  account  ol  the  I 1 uid-st ruclure 
interaction.  An  interpolation  sclieme  was  used  to  calculatt'  1 lu’  ci'itical 
loads  between  the  lon^-durat i on  and  short -dural  ion  extremes.  While  this 
theory  at^rced  reasonal)ly  well  with  t lie  test  results  on  wiiiih  the  Kxcess 
Impulse  Kule  is  based,  it  still  did  not  provide  a sa  t i s I ac  t orv  t heori-- 
lical  basis. 

A second  attempt  resulted  in  a theorv  that  dii'erued  substantially 
Irom  the  lest  ri’sulls  on  which  I lie  Kxcess  Impulse  Hiile  is  based.  At  this 
point  we  tried  lo  ( itid  more  details  on  I liosc  e x[ie  r i men  I s to  I rv  to  deler- 
mine  when'  the  theory  was  wronp:--in  the  representation  ol  the  structure 
or  in  the  model  ol  the  I 1 u i d -st  rue  t u re  interact  Ion.  We  were  unable'  lo 
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find  more  details  on  those  experiments  and  beuan  a seareli  lor  other  ex- 


I 


perimental  results,  particularly  scale  model  experiments. 

We  eventually  found  a few  model  tost  results  1 rom  1951-55  and  per- 
formed a limited  number  of  confirmatory  tests  in  our  latioratory.  Tliese 
results  showed  that  buckling  was  the  predominant  response  mode  ol  t he 
pressure  hull.  Hence  an  attempt  was  made  to  develop  a i-evisod  theory 
based  on  buckliiif>-,  instead  of  simple  rigid-plastic  response. 

The  work  presented  hero  combines  dynamic  buckling  ol  cylindi’ical 
shells,  developed  at  SRI  over  the  past  15  years,  with  a Russian  treat- 
ment of  fluid-structure  interaction.  The  result  is  a plausible  theory 
based  on  experimental  observations.  W’e  believe  it  [joints  the  way  to 
theoretical  understanding  of  otlior  rcs|3onse  modes  involving  buckling  and 
flu  id -St  rue t ure  interaction. 
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2.  THEORETICAL  APPROACH  .WD  liESULTS 

In  tlie  tlicorv  doveloped  here,  dynamic  overall  buckling  is  taken  as 
the  critical  mode  lor  hull  damage.  Experimental  results  show  that  the 
damage  gradient  with  increasing  load  is  steep;  lienee  an  analysis  based 
on  threshold  response  is  adequate. 

The  theory  describes  buckling  caused  by  a transverse  incident  pulse 
in  water  with  specilic  results  lor  a rectangular  pulse,  that  is.  a pulse 
with  an  instantaneous  rise,  a constant  pressure,  and  an  instantaneous 
decay.  The  rectangular  pulse  is  a reasonable  idealization  ol  a shock 
wave  from  an  underwater  nuclear  explosion,  with  surlacc  culoll. 

The  wavelengths  observed  in  experiments  indicate  that  the  predomi- 
nent  Inickling  mode  forms  during  the  initial  clastic  phase  of  de  f o niia  t i on  . 
This  mode  enters  the  initial  conditions  ol  the  plastic  phase  and  is  the 
prelcrrcd  mode  lor  lurther  amplification.  Therefore,  the  theoretical 
approach  is  to  treat  the  initial  response  of  the  shell  by  elastic  modal 
analvsis  until  plasticitv  occurs  in  the  hoop  mode,  and  to  use  the  dis- 
placement and  velocity  d i st  ri  Itu  t iiai  s at  this  time  as  initial  conditions 
for  t lie  subsequent  plastic  restionse.  .Amplification  ol  the  preferred 
tnickling  mode  is  assumed  to  continue  during  the  plastic  phase  until  the 
hoop  mode  reaches  the  maximum  inward  rarlial  d i sjil  acemen  t . The  experi- 
mental and  theoretical  final  delormcd  shapes  are  similar,  indicating 
the  validity  of  the  approach. 

It  is  also  observed  in  exiierlmcnts  that  a central  length  of  shell 
buckles  to  a shape  that  rloes  not  depend  on  the  axial  position  of  the  cross 
section.  This  rieforniat  ion  pattern  is  seen  in  the  shell  shown  in  figure  2.1. 
Also,  the  cross  section  of  the  shell  in  figure  A . 1 1 has  a shape  that  is 
tvpical  ol  the  shapes  of  slices  that  were  cut  along  the  shell  for  a 


cenlrnl  leiiRth  covoring  several  diameters.  This  observation  justified  the 
simplityins  assumption  of  analyzing  dynamic  buckling  of  a I'ing  instead 
ol  a cylindrical  shell  of  finite  length,  a much  more  complicated  problem. 

Experimental  and  theoretical  pressure  histories  at  the  lirst  con- 
tact ol  tlic  pulse  with  the  shell  also  show  good  agreement  exhibiting  a 
sliai’i)  rise  lollowed  by  a sharp  decay  and  gradtial  I'ise  (Figures  3.12  and 
5.1(1).  The  final  ilecay  dilfei’s  because  tlie  tlieory  is  based  on  an  ideal- 
ized rectangular  pulse.  Aside  from  flic  1 inal  decay,  the  main  difference 
in  the  pulses  is  the  pi'ossurc  level  in  the  valley  between  the  sliarp  de- 
cay and  gradual  rise.  This  pressure  level  depends  on  the  nature  ol  the 
structure,  and  since  the  structiires  differ,  tliese  levels  slioultl  not  be 
the  same . 

An  eiiective  may  oJ  j'epresejil  ing  tlie  datiiage  ol  a structure  to  jvtilse 
loailing  is  by  an  isodamage  curve  in  tlie  P-1  (pressure- impul  sc)  plane 
(Figure  6.2).  For  an  incident  rectangular  pulse  ol  pressure  P and  dura- 
tion T,  tlie  imimlse  is  I = PT . Damage  to  a submarine  shell  may  be  taken 
as  the  maximum  peniianent  radial  displacement  ; in  dimensionless  tenns 
this  displacement  may  be  represented  as  a fraction  of  the  shell  ratlins 
Ol'  a multiple  t)l  the  shell  tliickness.  A pressure- impul  se  curve  is  there- 
Itji'o  tile  curve  in  the  P-1  jilaiie  on  which  • is  a constant;  the  curve  con- 
nects rectangular  pulses  ol  etiual  el  lectiveiiess  that  range  1 rom  itleal 
impulse  to  step  loading. 

The  theory  developed  here  Was  used  to  generate  P-1  curves  Itir  an 
externally  stillened  cylindrical  shell.  The  pressurt'  and  duration  til 
the  incident  rectangular  pulses  are  represented  bv  P and  f , and  the 
damage  is  represented  nond i men s i on  a 1 I y by 

5 = max|w(0,t  ) - w (t  ) cost!  I 

a ' t II  ' 


- 


that  is,  l)v  the  maximum  inward  radial  displacement  at  the  time  t^  when 

the  hoop  mode  inward  displacement  w (t)  reaches  a maximiun.  The  displace- 

o 

ment  parameter  ^ is  the  linal  maximum  inward  displacement  divided  t)v 
a 

the  shell  radius  a. 

For  the  model  cylindrical  shell  subjected  to  a roctanftular  incident 

pulse  witli  a duration  efp.ial  to  a transit  time  ot  two  diameters  (t  ^ = 

100  sec,  T = 1),  an  incident  pressure  of  p.  = 2016  psi  (P  = 0.0062)  will 

cause  a peimianent  detormation  ot  2.5  wall  thicknesses  = 0.01);  the 

a 

corresponding:  incident  impulse  is  201  psi-msec  (I  = 0.0218).  The  static 

vield  pressure  is  p = c (k  H)/a  = 925  psi,  uivinn  a pressure  ratio 
y y a 

p /p  = 2.21;  thus  the  excess  impulse  is  I = (p  - p )t  = 112  psi-msec 
i y X t y d 

( I = 0.0136) . 

As  another  example,  if  the  pulse  has  a transit  time  ol  tour  diameters 

(t  = 200  Lisec,  T = 8)  , a pressure  of  p = 1250  psi  ( P = 0.0038)  will  also 
d i 

cause  a permanent  detormation  of  3.5  wall  tliickncsses  ('  = 0.01);  the 

a 

pressure  ratif)  is  now  p p = 1.35.  Tlie  incident  imttulse  is  250  psi- 

i y 

msec  (I  = 0.0301)  and  the  excess  impulse  is  I =65  psi-msec  (I  = 0.0075)). 

,x 

Fxamination  ol  the  isodamatte  curves  sliows  that  t lie  excess  impulse 
depends  on  the  P-I  or  P-'t  point  selected.  Hence,  excess  impulse  is  not 
constant  lor  constant  dnmay;e  and  is  tfiercfore  not  as  accurate  a measure 
of  the  dama>te  potential  of  underwater  shock  waves  as  llie  P-1  representation. 
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3.  EXPEKIMENTAL  BASIS  I OK  THEORY 


3 . 1 BACKGROUXU 


Model  tests  wore  conducted  to  simulate  the  response  modes  that 
occur  in  submarine  pressure  hulls  loaded  by  underwater  shock  waves  I rom 
large  explosions.  The  results  ol  those  tests  tormed  the  basis  tor  our 
theory . 

We  first  looked  tor  large-scale  test  results  to  determine  response 
modes.  However,  the  available  results  wore  very  sketchy,  undoubtedly 
because  the  damage  gradient  with  increasing  load  is  very  steep,  as  sliown 
below.  Because  of  this  steep  damage  gradient,  acliicvoment  o(  threshold 
response  requires  either  pi-ocise  critical  load  predictions  and  precise 
load  control  or  a largo  number  ot  tests.  Largo  c-xplosions  generally 
produce  loads  that  are  uncertain  by  at  least  10  percent,  wliich  is 
enough  to  span  the  i-ango  from  no  damage  to  extensive  damage.  Hence, 
except  for  an  occasional  luck  result,  many  test  structures  must  bo 
fielded  witti  large  explosions  to  got  threshold  response.  Since  this 
is  costly,  there  are  not  enough  large  scale  tests  to  illustrate 
threshold  response  modes. 


Conlrarv  to  our  expect  ions,  we  also  (ound  verv  few  modc'l  lest  re- 
sults to  review.  Intpiiries  to  the  Navy  led  to  I'n  i t ed  States  Na\al 

Research  atid  be ve lo|)mcn t Cetiter,  t'arderock,  as  1 lu'  best  source.  In  a 

2 , 3 

visit  to  ( arderock,  we  acttuiretl  two  referitices  dated  I'l.l  1 atid  ltt55, 
which  ilescribe  model  ti'sts  i ti  prt’partit  ion  lor  the  Pact  tic  nuclear  lest 
i\IfAV\M.  Xppari'iitlv,  no  other  model  tests  showing  threshold  ri'S|tonsi‘ 
have  l)een  reported.  fo  \eritv  tliese  results  and  obtain  (plan  1 i t a t i \ i- 
ri  spouse  inlormalion,  we  supp  I emeii  I ed  the  available  results  w i t li  some 
mode  1 tests  at  SRI. 
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Model  tost  results  illustrating  the  various  threshold  ros|)onse 
modes  evoked  by  simulated  underwater  shock  waves  from  nuclear  explosion 
are  shown  in  Figures  3.1  to  3.8.  The  response  modes  may  be  class! lied 
as  the  following  variations  of  hoop  response  and  a.xial  response; 

Hoop  response 

Quasi-static  overall  buckling  (Figures  3.1  and  3.2) 

Dynamic  plating  buckling  (Figure  3.3) 

Dynamic  overall  buckling  (Figures  3.1  and  3.5) 

.Axial  response 

Dynamic  plating  buckling  due  to  axial  load  (Figures  3.6 
and  3.7) 

Dynamic  plating  buckling  due  to  overall  bending  (whipping) 
(Figure  3.8). 

Mull  dimensions  arc  given  in  Figure  3.9  for  Figures  3.1  and  3.5  through 
3.8,  and  in  Figure  3.10  for  Figures  3.2  through  3.1.  Kach  of  tliese 
modes  and  the  mechanisms  by  wliich  they  are  activated  arc  discussed 
be  1 ow  . 


3.2  HOOP  RESPONSE 


3.2.1  Quasi -St  at ic  Overall  Buckling 

This  response  mode,  shown  in  Figure  3.1,  consists  ol  a single 
inward  lobe  extending  over  several  stiffeners.  In  this  example,  the 
lobe  is  about  two  hull  diameters  long.  .A  similar  model  loaded  hydro- 
statically is  shown  in  Figure  3.2.  The  response  is  essentially  the 
same  as  in  the  model  loaded  dynamically.  Hence,  the  term  quasi-static 
overall  buckling.  This  is  the  well-studied  classic  buckling  bcliavior 
ol  a finite  cylinder  under  external  pressure.  Buckling  occurs  wlicn  the 
hoop  thrust  in  the  hull  becomes  great  enough  to  overcome  the  hoop  bonding 
resistance  in  the  critical  mode.  The  nonun i formi 1 v of  the  dvnamic  load 
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docs  not  affect  the  basic  response  mode,  but  probably  affects  the  load 
recpiired  to  produce  a given  response  level. 

Quasi-static  overall  buckling  is  a special  case  of  dynamic  overall 
buckling  described  below. 

3.2.2  Dynamic  Plating  Buckling 

This  response  mode,  shown  in  figure  3.3,  consists  ol  inward  and 
outward  dimpling  of  the  fulll  plating  between  stiffeners.  Inward  dimpling 
occurs  alone,  but  outward  dimpling  always  occurs  between  two  inward 
dimples.  Dimpling  occurs  when  t fie  lioop  ttirust  is  sufficient  to  overcome 
tlic  tiending  resistance  ol  the  plating  in  ttiis  mode.  Since  this  response 
mode  results  in  relatively  minor  hull  damage,  it  is  not  considered  the 
critical  damage  mode. 

3.2.3  Dynamic  Overall  Buckling 

This  response  mode  i.s  shown  in  f'igui’cs  3.1  and  3.5.  It  is  simi- 
lar to  the  quasi-static  overall  tiuckling  in  tliat  the  lolios  extend  over 
several  stiffeners,  fnit  it  differs  in  ttiat  tfiere  are  several  lottos  in- 
stead of  one  (particularly  evident  in  Figure  3.1).  The  sliift  from  a 
single  lofie  to  several  lobes  in  flic  thrcstiold  mode  is  due  to  t fie  dil- 
lerenco  in  load  (iuration.  The  response  time  lor  t tie  development  of  a 
single  lotie  is  greater  t fian  for  several  lobes.  Hence,  file  load  duration 
mav  tie  sfiortor  ttian  necessary  evoke  the  single  lofic  mode,  fiut  long 
enough  to  evoke  higfier  order  modes.  Tfiis  type  of  dynamic  tiuckling  in 
higtier  order  modes  is  callect  pulse  tiuckling  and  has  lieen  studied  exten- 
sively at  Stfl  over  the  past  15  years  (see  Section  1 lor  references). 

Ttie  analysis  yields  lime-dependent  amplification  factors  lor  t tie  various 
modes.  Tvpicallv,  t tie  ampl  i I i cat  inn  factors  stiow  a peak  ttiat  extends 
over  a few  modes.  If  file  load  duration  is  long  enougii , t lie  peak  occurs 
in  file  region  ol  the  static  tiuckling  mode.  The  most -magn  1 1 it'd  modes 
predictttl  tiv  t lie  ttieorv  generally  agree  well  witli  experimental  results. 

if) 


3.3  AXIAL  RESPONSE 

3.3.1  Dynamic  Plating  Buckling;  Due  to  .Axial  Load 

This  response  mode,  shown  in  Figures  3.6  and  3.7,  consists  ol  an 
inward  deformation  of  the  plating  between  stiffeners  that  may  extend  all 
around  the  circumference  (Figure  3.6)  or  only  a few  stilfener  spacings 
( F igurc  3.7). 

3.3.2  Dynamic  Plating  Duckling  Duo  to  Overall  Bending  (Whipping) 

This  response  mode,  shown  in  Figure  3.8,  consists  of  inward 
buckling  of  the  plating  between  the  central  stiffeners  due  to  axial 
bending.  The  buckles  occur  on  the  compressive  side  and  extend  over  a 
few  stiffener  spacings. 

3.J  D I. SCISSION 

The  response  modes  described  abo\c  involve  instability  under  [julsc 
loads.  They  can  all  bo  treated  by  pulse  t)uckling  tlioory.  The  tlioory 
developed  hero  is  restricted  to  the  hoop  response  mode  ol  dynamic 
overall  tnickling,  which  is  taken  as  the  critical  hoop  mode  tor  hull 
damage.  As  shown  below,  the  damage  gradient  with  increasing  load  is 
steel);  lienee  an  analysis  based  on  threshold  response  is  adeciuate.  The 
theory  lor  dynamic  overall  buckling  developed  earlier  lor  loads  that 
are  independent  ol  the  response  was  e.xtondod  to  account  lor  tlie  effects 
ol  I luid-st ructure  interaction. 

The  general  goal  in  tiie  design  ot  stilleiicd  shells  lor  hulls  is 
to  match  the  critical  load  for  the  plating  lictwecn  st it  loners  and  the 
critical  load  lor  overall  bending  ol  the  plating  and  stilleners  together. 
Moi'eover,  t lie  stilleners  arc  designed  not  to  lail  by  lateral  buckling  as 
shown  in  figures  .i.d(b)  anrl  3.  1(a).  Hence  an  ideally  designcrl  stiflencr 
shell  woulfl  behave  essentially  as  a unilorm  shell  with  increased  bending 
St  i I fness . 
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To  examine  the  response  ol  unilorm  shells  to  underwater  shock  wave 
loading,  we  tested  a scries  ol  lonj;  alimiinum  (6061-T6)  cylinders  ol 
3-inch  diameter  and  1 16-inch  wall.  Cross  sections  ol  two  ol  the  cylin- 
ders are  shown  in  I'iKure  3.11.  At  liif^her  loads  the  central  portion  on 
the  loaded  side  is  di’ivcn  inwarti,  resulting  in  rupture  alon^  the  loaded 
side.  The  cross  sections  ol  li^uro  3.11  sliow  that  the  mode  numbers  ol 
interest  in  the  threshold  region  are  ol  the  order  ol  oiRht.  This  observa- 
tion was  used  in  dcvelopinj;  the  theory. 

Figure  3.12  shows  the  pressure  history  mcasuretl  on  the  surlace  oi 

a model  like  those  ol  Figure  3.11.  This  record  was  ol)tained  witli  a 

thin  toil  (~  1 mil)  ytterbium  f;afic  attaclicd  to  the  surlace  ol  I lie  model 

•1 

at  the  point  ol  lii'st  loading;.  .\  similar  record  obtained  by  other 
moans  is  shown  in  .\ppendix  ,\  (Fi};iire  ,\-3). 

II  dimpling  ol  the  liull  plating  (Fif;iirc  3.3)  occurs  !>clore  ovei’all 
bucklinfi,  it  will  inilucnce  the  overall  l)ucklinn  mode  l)y  introducing 
nonun  i lormi  t i cs . Thus,  tlie  tlicory  described  in  Section  1 should  bo  re- 
f^arded  as  a lirst  sto|)  in  applying  pulse  bucklinp;  theory  to  .•submarine 
pressure  hulls,  takin>;  account  ol  the  ellects  >1  t 1 uid-st  ruct  urc- 
in  teract ion . 

•As  mentioned  at  the  bcf^iTinint;  of  tliis  sectic'n,  t lie  damatte  itradicnt 

in  underwater  structural  res|)onse  is  steep.  The  liasis  tor  tliis  oliserva- 

I ion  is  fiiven  in  Figure  3.13.  llie  ordinate  in  1 inure  3.13  is  1 lu'  ratio 

of  tile  peak  del  lection  A to  the  wall  thickness  t . The  aliscissa  is  t lie 

i-atio  ol  the  load  ( livd  ros  i a t i c pressure’  P plus  (leak  I rci‘  liehl  shock 

o 

pressure  P ) to  t tic  colIa|ise  pressure  ol  the  liull  P . U'e  see  tliat  t lu' 
m c 


del  lection  increases  from  O to  more  t tian  10  wall  t hi  cknessis  tor  a 10 
percent  increase  in  load. 


Comparison  of  figures  3.3(a)  and  (b)  with  3.4  further  illustrates  the 
steep  damage  gradient  for  underwater  structural  response.  Hero  an 
increase  in  the  free-field  pressure  from  750  to  1050  psi  produced 
response  ranging  from  threshold  to  hull  rupture. 

The  steep  damage  gradient  for  underwater  structural  lesponse  means 
that  the  theory  for  critical  load  curves  can  be  based  on  threshold 
rosi)onse.  This  is  very  important  because,  in  the  analysis  of  threshold 
response,  we  can  take  advantage  of  many  of  the  simplifications  and 
idealizations  of  classical  mechanics. 
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FIGURE  10b  - MODEL  LH-7  AFTER  TEST  3-C 

THE  FAILURE  IS  OF  THE  OVERALL  INSTABILITY 
TYPE 


MP  1714  61 

FIGURE  3.1  QUASI-STATIC  OVERALL  BUCKLING  PRODUCED 
BY  AN  UNDERWATER  EXPLOSION  (FIGURE  10b 
OF  REFERENCE  2) 


Ml'  1 7)4  02 


FIGURF  3 2 QUASI  STATIC  OVERALL  BUCKLING  PRODUCED 
BY  A HYDROSTATIC  LOAD  (FIGURE  1?  OF 
REFERENCE  31 
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IrfI  DYNAMIC  PLATING  BUCKLING  PRODUCFD  BY  AN  UNDERWATER  EXPLOSION 
AND  BACK  VIEWS 


dynamic  plating  buckling  produced  by  an  underwater 

EXPLOSION  TOP  AND  BOTTOM  VIEWS 


MOUHE  3 3 


CHARGE 


I 

DETONATOR 


PRESSURE  TRANSDUCER  LOCATIONS 


CALIBRATED  PEAK  PRESSURE  VERSUS  STANDOFF 


Idl  RESULTS  OF  CHARGE  CALIBRATION  TESTS.  MODEL  OF  FIGURES  3 3a  AND  b 
WAS  AT  24  INCH  RANGE.  MODEL  OF  FIGURE  3.4  WAS  AT  18  INCH  RANGE 

VA  1714  3A 


FIGURE  3 3 (Contituied) 


PRESSURE  transducer  4 


PRESSURE  TRANSDUCER  2 


PRESSURE  TRANSDUCER  3 


l«l  TYPICAL  PRESSURE  RECORDS  FROM  CHARGE  CALIBRATION  TESTS 
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FIGURE  lOf  MODEL  U-b  AFTER  TEST  2 C 


FIGURE  10c  MODEL  LH  9 AFTER 
TEST  1 C 


THE  INTERNAL  LEAD  KEEL  IS  VISIBLE  THE 
KEEL  AND  PART  OF  THE  PLATING  WERE 
ROTATED  ABOUT  ONE  EDGE  OF  THE  KEEL 
IN  THIS  TEST. 


THE  TOTAL  PEAK  PRESSURE  WAS 
ABOUT  THE  SAME  AS  FOR  LH  3 
IN  FIGURE  lOd  BUT  THE  DEPTH 
WAS  GREATER  BY  314  FT 


FIGURE  3.5  DYNAMIC  OVERALL  BUCKLING  (FIGURES  10. 
REFERENCE  ?) 


FIGURE  3.6  DYNAMIC  PLATING  BUCKLING  DUE  TO  AXIAL  LOAD 
(FIGURE  10a  OF  REFERENCE  2) 


(b)  FIGURE  lOe  MODEL  A 1 AFTER  TEST  2 C 

THE  DISHING  AT  THE  ENDS  IS  MORE 
EXAGGERATED  THAN  IN  FIGURE  lOa 


FIGURE  3.7  DYNAMIC  PLATING  BUCKLING  DUE  TO  AXIAL  LOAD 
(FIGURE  lOf  OF  REFERENCE  21 


(cl  FIGURE  lOd  MODEL  SH  7 

THE  CENTRAL  BAY  SHOWS  FAILURE 
ASSOCIATED  WITH  SHIRRING  THERE  ARE 
LOBES  IN  ADJACENT  BAYS 


FIGURE  38  DYNAMIC  PLATING  BUCKLING  DUE  TO  OVERALL  BENDING 
(WHIPPING)  (FIGURE  lOd  OF  REFERENCE  2) 
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FIGURE  3.9  PRINCIPAL  DIMENSIONS  OF  MODELS  OF  FIGURES  3?  AND  3.5-3  8 (FIGURE  1 OF  REFERENCE  2) 


FIGURE  311  DYNAMIC  OVERALL  BUCKLING  OF  UNIFORM  SHELLS  (Al  6061 -T6,  3 INCH 
DIAMETER,  1,16-INCH  WALL)  PRODUCED  BY  UNDERWATER  EXPLOSIONS 
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FIGURE  3 12  SURFACE  PRESSURE  AT  POINT  OF  FIRST  LOADING  ON  A MODEL  LIKE 
THOSE  OF  FIGURE  3.11,  AS  MEASURED  WITH  A THIN  FOIL  (-1  rnil) 
YTTERBIUM  GAGE  ATTACHED  TO  THE  MODEL  (REFERENCE  4) 
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FIGURE  3.13  maximum  PERMANENT  DEFORMATION  VERSUS  RATIO  OF  TOTAL  PEAK 
PRESSURE  TO  HYDROSTATIC  COLLAPSE  PRESSURE  FOR  MODELS  OF 
SUBMARINE  PRESSURE  HULLS  (FIGURE  9 OF  REFERENCE  21 


■1. 


BUCKLING  THEOHV  OU  SL15MERGKD  CVLINDIilCAL  SHELLS 


•1 . 1 BACKGROUND 

Because  cat  ast  roplii c daniafie  occurs  when  the  cylindrical  shoot  and 
ring  stitfonors  buckle  together,  as  shown  in  Figure  3.1(a)  ol  Section  3, 
the  theory  was  developed  to  doscriljo  overall  buckling  rather  than  the 
local  ttuckling  that  consists  only  ol  sheet  dimpling,  shown  in  Figure 
3.3(a).  The  theory  describes  buckling  caused  by  a transverse  incident 
pulse  in  water  with  s[>ocific  results  lor  a rectangular  pulse  (instantane- 
ous rise,  a constant  pressure,  and  an  instantaneous  decay).  The  rectan- 
gular pulse  is  a reasonable  idealization  of  a shock  wave  from  a large 
underwater  explosion,  with  surface  cutoff. 

The  wavelengths  observed  in  the  experiments  indicate  that  the  pre- 
dominant buckling  mode  forms  during  the  initial  elastic  phase  of  deforma- 
tion; thus  this  mo<lc  enters  the  initial  conditions  of  the  plastic  phase 
and  is  consequently  the  prcferrefl  mode  lor  further  am|jli  I icat  ion  . Tliere- 
fore,  the  theoretical  approach  is  to  treat  the  initial  response  of  the 
shell  by  elastic  modal  analysis  until  plasticity  occurs  in  the  hoop  mode, 
and  to  use  the  displacement  and  velocity  distributions  at  this  time  as 
initial  conditions  lor  the  subsequent  plastic  response.  .Amplification 
of  the  prelerred  buckling  mode  is  assumed  to  continue  during  the  plastic 
phase  until  the  membrane  mode  reaches  the  maximum  inward  radial  displacement 

The  buckling  analysis  is  based  on  the  dynamic  pulse  buckling  theory 

developed  at  SRI  over  the  past  15  years.  The  f luid-st  nictu re  interaction 

5 

is  based  on  the  Russian  work  by  Mnev  and  Pertsev.  The  relevant  parts  of 
this  work  are  included  below  to  make  the  treatment  complete.  .Appendix  B 
summarizes  the  theory  presented  in  this  section  in  the  form  used  lor  I lie 
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computations  reporiod  in  Sections  5 and  6.  Appendix  B also  includes 
several  impro\ed  approximations. 


1.2  GOVEIIMNG  EQUATION  OE  INITIAL  ELASTIC  PHASE 

It  is  assumed  that  tlie  cylindrical  shell  is  long  enough  for  a 
central  portion  to  respond  as  an  infinitely  long  shell.  The  theory  is 
thcrelore  planar  and  is  described  in  polar  coordinates  r,0. 

Considering  unit  axial  length  regarded  as  a ring,  tlie  proportionality 
of  moment  and  curvature  cliangc  is 


M = Elx 


(1.1) 


where  E is  the  clastic  modulus  of  the  shell  material,  I is  tlie  second 
moment  of  area  of  the  ring  cross  section,  and  x is  the  curvature  increase 
corresponding  to  a radially  inward  displacement  w(H,t).  The  expression 
f o r X is 


1 / 0 w 


(1.2) 


where  a is  the  radius  of  Uic  centroid  ot  the  ring  cross  sectional 


Figure  1,1  shows  a shell  element  at  time  t that  has  l)cen  dis- 
placed radially  inwards  a distance  w(0,t),  shortened  from  an  original 
lengtti  of  adU  to  dl,  and  fjont  so  ttiat  t fie  sid)tcndcd  angle  clianged  from 
clH  to  d . Tliese  del  oniiations  were  caused  lj\  the  external  pressure  p(B,t) 
and  were  resisted  t)y  f he  internal  mcmtirane  thrust  .S , bending  moment  M, 
and  sliear  force  t),  along  wifh  I tie  radial  inertial  force.  .As  an  apnroxi- 
mation,  the  ci  rc'umlerent  ial  displacements  and  variations  in  S are 
neglected.  Neglecting  also  the  rotary  inertia,  the  eriuations  <il  motion 


n ( t he  e lemc n t a re 
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(1.3) 


(1.1) 


whore  m is  the  mass  per  unit  arc  length  ol  shell.  In  (1.3),  the 
curvature  is 


= - ^ X (1.5) 

3'  a 


.-\  further  approximation  is  made  by  sotting  d).  = adO  . 

Equations  v-1.1)  through  (1.5)  are  now  combined  to  produce  the 
governing  equation  in  the  form 


E I / 3 w 3 w^ 

■1  \ 1 ^ 
a \3B  3U'7 


1 / 3 w 

2 \ 2 
a V30 
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3 w 
•? 

3t“ 


(1.6) 


The  thrust  S is  the  product  of  the  av^erage  c i reunite  ren  t i a 1 com- 
pressive stress  c and  the  cross-sectional  area  .A  per  unit  length  ot 
shell.  If  H is  the  depth  of  the  stiffener,  ns  shown  in  the  two  cases 
in  Kigurc  1.2,  the  area  .A  may  be  expressed  as 


A = k H 
a 


(;  * t - n^' 


(1  .7) 


where  ' - II  h and  \ . 1,  ft  . .As  indien'ed  in  Figure  1.2,  I,  is  the 

stiflener  spacing,  I is  the  stiffener  thickness,  and  fi  is  the  slu’Ct 
thickness.  Bec..','se  the  membrane  stress  is  assumed  to  be  approx  ima  t e 1 v 
constant  around  the  e i rcum  1 1>  rence  at  ettch  instant,  it  tolt>.iws  that 
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Ew  a,  where  E is  the  elastic  modulus,  e 
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is  the  avcrat;e 
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c i reunite  rent  ial  compressive  membrane  strain,  and  w 


radially  inward  displacement.  Thus 


S = k HEw  a 
a o 


( l.K) 


For  the  cross  sections  ot  Figure  T.2,  the  second  moment  ol  area  may  bo 
expressed  as 
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I = k H , 12 
i 


C ^ i - 1)  - (i  - DC  - 1) 


/\:  c ^ - 1) 


(1.9) 


It  the  density  of  the  shell  material  is  p , the  mass  per  unit  length  is 

s 


m = 0 A = p k 11 

s s a 


(1  . 10) 


Uy  substitution  of  S,  I.  and  ni  from  (4,8),  (4.9),  and  (1.10)  into 
(1.6),  the  Kovernitift  ecpiation  becomes 
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Note  that  the  coefficient  w in  this  partial  differential  equation  is 

o 


a function  of  time,  as  is  the  pressure  p . 

•As  an  alfjebraic  simpli  f icat  icn , dimensionless  quantities  arc  intro- 
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to  Table  1 

1.1,  to  put 
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The  t inal  1 Orm  ol  the  fiovcrnin^  cM|uat  ion  can  be  oljtaineil  onlv  by 
takinti  into  account  the  depettdence  ot  the  pressure  on  the  response  ol 
the  shell.  This  1 lu id-st rucl ure  interaction  is  the  subiect  ol  the  next 
sect  ion . 


•1.3  ILUID-SHKLL  INTEiaCTlON 

5 

In  the  Work  ol  Mnev  and  F’ertscv,  the  pressvire  I ield  in  the  Iluid 
is  considered  as  the  superposi  t i<tn  of  three  1 iolds:  (1)  pressure  in  the 

undisturbed  incident  pulse,  (2)  pressure  caused  by  ditlractif)n  of  a I i.xed 
riftid  cylinder,  and  (3)  pressure  caused  by  radiation  from  a deforminfi  and 
translating  cylinder.  This  superposition  is  justified  in  the  tollowinR 
0 xpos i t ion . 


Let  ^ be  the  velocity  potential  lor  the  diffraction  field.  Then, 
the  par  icle  velocity,  pressure,  and  wave  ecpiat ion  are. 


V 

d 


a 

d 


3r 


(1.13) 


bet  V be  the  radial  comnonent  ol  ve loci  tv  on  a fixed  unit  circle 
s 

(r  = 1)  associated  with  the  passintt  undisturbed  incident  wave.  Foi-  a 
fixed  ri^ii<l  cylinder  to  occupy  the  unit  circle,  the  incident  and  dif- 
Iraction  fields  musi  superimpose  to  provide  no  radial  velocity  on  the 
unit  circle';  that  is. 


(1.11) 


The  solution  ol  I lu*  cylindrical  wave  '.tpiation  in  ( 1.13),  satislyinu  the 

boiindarv  condition  (1.11)  on  the  unit  circle  and  the  condition  (I 

d 

as  r ♦ » (initial  conditions  = (I,  9-  9t  = (i)  mav  lie  found  by  a 

d d 

Laplace'  transformation  with  re'spe'ct  to  time  t and  separation  ol  tarialile's 
with  ri'spLe-t  to  the'  polar  coordinate's  i',(). 
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It  the  radial  component  ot  the  velocity  in  the  incident  wave  at  the 
unit  cii’Cle  has  the  Fourier  representation 


■ = Z 

s .t— ( n 


(t)  cosnO 


the  Laplace  transform  ot  the  ditfraction  field  velocity  potential  is 


I)  K ( rs ) . cf)snH 

Ln  n 

sK  (s) 

n=0  ^ n-1 

K ( s ) n ; — ; — 

n k ( s ) 

n 


In  ( 1.1(5),  s is  the  transtonn  parameter  and  K (s)  is  a moditied 

n 

Bessel  function  ot  the  second  kind  of  oriler  n and  argument  s.  Trans- 

L L 

tormint;  the  prcssui'e  tomuila  in  (1.13)  t;ives  p = -s  - and  hence,  l)v 

d d 


(1.  1(5)  , 


' sb^  m cosnO 
n n 


wliere 


sK  ( s) 
n-1 

K (s) 
n 


K (s) 
n 

sK'  ( s ) 
n 


If  (lie  inciefent  wave  pressure  evaluaterl  on  t fie  unit  circle  i .s 


p = Z 

s ^ n 


and  the  pressuro  distribution  on  the  *ixotl,  rit;id  cylinder  is 


o o 

= 2- 


(1.20) 


II 


it  follows  that 


oL 

P 


^ (a  - sb  m ) c<jsnH 
s ci  ^ ' n n n / 

n = 0 


(1.21) 


and 


oL  L L 

p = a - sb  m 
n n n n 


(1.22) 


Inversion  of  (1.22)  provides  the  exact  solution  for  the  pressure  on  a 
fixed,  ri^id  cylinder  suljjectod  to  a transverse  pressure  pulse,  but  it 
is  more  fruitful  to  employ  the  approximate  method  to  be  described. 

Let  ^ be  the  velocity  potential  lor  the  radiation  field.  Then, 


9" 
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V r 
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■■ 

Recallins  that  the  radially  inward  displacement  ol  the  dolormable  shell 

boundary  of  the  fluid  is  w(0,t),  the  benmdary  condition  to  be  satisfied 

by  2 is 
r 


9w 


0 


r - 1 


( 1 


.2  1) 


II  the  shell  d 1 S|)l acemen t is  niven  the  Tourier  representation 


act 

w = 7 w cosnH  (1.25) 

< n 

n n 


the  t ran s 1 ormed  bound  a rv  condition  (1.21),  with  the  initial  condition 
w(H,n)  = (),  is 


L 

sw  cosnO 


( 1 . 25  ) 


The  solution  of  the  cylindrical  wave  equation  of  (1.23)  satisfying  the 


boundary  condition  (1.24)  on  the  unit  circle  and  the  condition  -►  0 

r 

as  r -♦  w (initial  conditions  = 0,  9"'  /9t  = 0)  may  be  found  by  a 

r r 

Laplace  transformation  with  respect  to  time  t and  separation  of  vari- 


ables with  respect  to  the  polar  coordinates  r,0.  This  procedure  leads 
to  the  transformed  radiation  velocity  potential 


00  sw^  K (rs)  cosnO 

n n 


sK  (s) 
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K ( s)  n 
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K ( s) 
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I'L.  L 

Transforming  the  pressure  formula  in  (-1.23)  gives  p = hence, 

by  (4.26) 


2 L 

= - / s w m cosj 

n n 


representat ion 


p = - 2.  p. 


it  follows  that 


P*"^  cosnH 


and , by  (4.27), 


rL  2 L 

p = s w m 
n n n 


Introducins  the  function 


!p  (S) 

n 


clef  inocl 


L 

Ip  (s)  = m ( s ) 

n n 


(1.31) 


e.xprcssion  (1,30)  becomes 


rL 
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2 L,L 
s V ip 
n n 


(4.32) 


With  initial  conditions  w (0)  = 0 and  w (0)  = 0,  whore  a dot  denotes 

n n 

time  differentiation,  the  transform  of  each  shell  acceleration  component 


is 


00 

2 L f ..  , -St 

sw  = lw(t)e  dt 
n J n 

o 


so  that,  ijy  applying;  the  convolution  theorem  of  transforms,  (1.32)  can  bo 


inverted  to  give  the  radiation  pressure  components 


t 

p*  = I w (t)^  ( t - r)dT  (1.33) 

n J n n 
o 


A physical  interpretation  can  be  given  to  the  function  ip  if  we 

th 

consider  the  unit  step  velocity  in  the  n mode  of  deflection, 

w (t)  = H(t),  whore  H(t)  is  the  unit  stop  function  (H  = 0 for  t < 0 

^ L -2 

and  11  =.  1 lor  t > (1)  . Thus,  w (t)  = tll(t)  and  w (s)  = s , which,  by 

n n 

t*Ij  ij  I* 

(I  32)  gives  t)  " = Ip  . Thus  = p is  the  radiation  pressure  asso- 

n ^ J n n n 

dated  with  the  n mode  for  a shell  moving  at  unit  velocity.  The 


f unct i on 


is  called  the  unit  function  of  resistance.  Other  proper- 


n 

ties  of  if  (t)  arc  ip  (0)  = 1 , («>)  = 0,  and 

n n n 
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(t)dt  = - 

1 n 


This  last  result  can  be  shown  to  be  the  apparent  addition  mass  during 
motion  of  a shell  in  an  incompressible  fluid. 

The  unit  function  of  resistance  can  also  bo  introduced  into  (1.22) 
for  the  transform  of  the  diffraction  pressure  components  so  that 


oL  L L L 

p = a - sb  ih 
n n n n 


(4.35) 


Inversion  of  (-1.35)  gives 


p = a ( t ) - b 
n n ] 


(t)  - \ b (T)ii 

n J n ^n 


(t  - -)dT 


(1.36) 


The  first  step  in  the  solution  of  equation  (1.12)  is  to  set 


o 1- 

p = Pp  + p 


(1.37) 


where  P is  a number  that  determines  the  pressure  in  the  incident  pulse. 

lor  P = 1 the  pressure  et|uals  the  bulk  modulus  of  water,  so  the  values 

-3 

of  P that  are  of  interest  arc  0(10  ).  The  second  stop  is  to  give 

o r 

w,  p , and  p their  lourier  scries  representations  (1.25),  (1.20),  and 

(1.23).  Kquating  coefficients  of  cosnH  >hen  gives  an  infinite  set  of 

ordinary  differential  equations.  The  third  step  is  to  approximate  the 

functions  ip  in  tlie  pressure  coefficients  (1.33)  and  (1.36)  anti  to 
n 

ev'aluate  the  pressure  and  velocity  coefficients  a and  b for  the 

n n 

given  incident  pulse. 
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■1.1  APPROXIMATE  UNIT  RESISTANCE  FL'NCTIONS 


Approximate  unit  resistance  functions  are  needed  because  the  exact 
functions  are  inconvenient  to  use  in  analysis.  The  choice  ol  the  approxi- 
mate functions  ip  ts  based  on  knowledge  of  the  exact  functions  (obtained 
n 

numerically  for  the  lower  modes),  general  properties,  and  the  need  lor  an 
analytically  tractable  form.  The  approximate  functions  chosen  to  repre- 
sent the  exact  functions 


♦„(.) 


= (1  - a) 


a = 0.S5,  t S 8 


(4.38) 


ip  (t) 
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('  -I'lh 


(4.39) 


These  arc  compared  with  the  exact  functions  for  n = 0,  1,  and  2 in 
Figure  1.3.  .An  improved  version  of  tp^  is  also  shown  and  is  used  in 
■Appendix  B. 

The  functions  (1.39)  are  approximations  in  the  form  of  two  straight 

linos  and  satisfv  the  conditions  ip  (0)  = 0,  Ip  (®)  = 0,  and  (1.34).  The 

n n 

function  (1.38)  is  special  because  two  straiglit  lines  do  not  form  an 

adequate  approximation  for  the  important  hoop  mode;  Ip^^(O)  = 1 but 

lp^^{=°)  = 1 - a,  so  a limiting  time  of  applicability  of  t ■^  8 is  suggested; 

when  t > 8,  set  ip  =0. 

0 

Substitution  of  ip^^  i”  (4.36)  and  (1.33)  gives  the  hoop  mode  pres- 
su  re  coe I f i c i on  t s 


p = a - b 
(1  o 
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(1.41) 
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i. 


and  substitution  oi  <p  in  (-1.36)  and  (1.33;  fjivos  the  other  pressure 

11 

coe  f f icien t s . 
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Improved  tomuilas  lor  p^^  and  p^  are  presented  in  .Appendix  B. 


■1.5  RKPRKSK.N'TATION  Of  A RECTANGUL.-\I!  PULSE 

With  the  aid  ol  UiKure  1.1,  it  is  seen  that  a rectam^ular  incident 
pulse  may  be  represented-  by  the  formula 


p = H(t  - 1 . cosO)  - H(t  - 1 ^ cosH)  (1.15) 

s 


the  time  orinin  bcinp  when  the  wave  front  first  contacts  the  unit 
circle.  In  (1.'15)  the  first  term  represents  step  loading  and  the 
second  term  represents  step  unloadinp:  after  a duration  T . The 
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corresponding  particle  velocity  I’esolv'ed  along  the  outward  normal  of 
the  unit  circle  is 


V = - cosB'HCt  - 1 + cosfi)  - H(t  - T - 1 + cosO)  ] (1.-46) 

s 


In  the  series  representations  (1.19)  and  (4.15)  for  p^  and  v the 
coefficients  are 


a = — I p dO 
o n 1 s 


a = \ p cosn0  dft 
n 1 s 


n a 1 


(4.47) 


b = - \ V dB 
o “is 


b = — 1 V cosnO  d0  n > 1 (4.48) 

n I s 


Substituting  (4.45)  and  (4.46)  for  p and  v in  the  integrals  (4.47)  and 

s s 

(4.48)  and  carrying  out  the  integration  give 


a = a'"H(t)  - H(t  - 2)'  ^ H(t  - 2) 


+ a "H(t  - T)  -H(t  - T - 2)"'  -H(t  - T - 2) 
n 


a "H(t)  - II(t  - 2) 
n 


a lUt 
n 


H(t  - T - 2)‘ 


b = b'^IKt)  - il(t  - 2)-  + b"[H(t  - T)  - H(t  - T - 2)] 
o O ■ o 


b = b' 'H(t)  - H(t  - 2)’  - H(t  - 2)  . b''''l!(t  - T) 
11  ' 1 


- H(t  - T - 2)'  . II(t  - T - 2) 


(4.49) 


b = b'"H(t)  - H(t  - 2)‘  * b"^H(t  - T)  - H(t  - T - 2); 
n n ■ n 
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where 


a = 
O 


a = — sin  nfi 
n nn  0 


n > 1 


— sinfl 


1 1 
- (fi  + - sin20  ) 
-02  0 


1 1 

sin(n  . 1)0  i sin(n  - 1)0 

n ^ 1 0 n - 1 0 


a = 
0 


a = - — sin  ti0 
n -n  1 


n > 1 


" 1 ^ 
b = — sin0 
0 - 1 


H \ 1 

h =— (0  sin20) 

1-12  1 


" 1 
b = - 
n 


1 1 

sin(n  - 1)0  - sin(n  - 1)0 

n . 1 1 n - 1 1 


COS0  =1-1 
0 


cos0^  = 1 - (t  - T) 


As  shown  in  riirure  1.1,  0 and  0 arc  the  an>iles  subtended  bv  the  wave 

()  1 

front  and  wave  rear. 


1.0  .SYSTEM  OF  CKIVEHNING  EQUATIONS  OF  ELASTIC  PHASE 

•After  substituting  the  series  representations  (1.25),  (1.20),  and 
o r 

(1.28)  for  w,  p . and  p in  the  governing  efpiation  (1.12)  and  equating 
to  zero  the  coefficients  of  each  cosnO,  wo  obtain  the  following  system 
of  equations; 

w 4 k w = a.(Pp  - p ) (1.50) 

o 2 o ' ()  O' 

2 2 or, 

w 4^k(n  - l)(kn  - w)w  = (Pp  - p I n '■  1 (1.51) 


1 


To  obtain  (-1.50)  it  is  assumed  that  w ■«  1 (hoop  displacement  small 

o 

compared  with  shell  radius).  Introducing  in  (1.50)  and  (1.51)  the  radia-  j 

tion  pressure  coefficients  (l.ll)  and  (1.13)  results  in  the  equations 


w + ..:w  + (k  - u,/2)w  = uP|J°  ( 1.52) 

o o 2 o 0 


(1 


:v  )w 
n n 


n n 


(n 


l)(k^n 


w ) 
o 


.O' 


w = (..Pp 
n n 


n > 1 (1.53) 


To  obtain  (1.52),  wo  noftloct  the  last  term  of  (l.  ll)  (this  term  is  re- 
tained in  the  siumnary  of  .\ppondix  B)  . The  accompanying  initial  condi- 
tions arc 

w (0)  = 0 w (0)  =0  n S'  O ( 1.5  1) 

n n 


■1.7  SOUTION  OF  GOVKUNING  EQUATIONS  OF  EL.-\STIC  PIL\SE 


The  solution  of  (1.52)  for  the  hoop  mode  with  initial  conditions 
(■1.51)  is 


w 

o 


t 


o 


So(  I -"T  ) 
C “ 


) 


o 

p„(-) 


d- 


(1.55) 


where 

= (v  - u.)/2 
s^  =-(v  ♦ l.)/2 

o = - Uk,^  - -.1/2)  ]“ 


£2i 


.53 


r 


i 


The  lirst  equation  of  1 lie  set  (T.53)  y;overnin^  shell  translation 


(1.56) 


and  lomulas  (-1.44)  give  the  coefficients 


a 


^4  = ' 


V = 0 
1 


f oi’  t ^ 2 


and 


O'  = 0 


= 0 


V = 1 
1 


tor  t > 2 


Hence  equation  (4.56)  is  replaced  by  the  followins  two  equations: 


w + uw  - (u/2)w  = '.^Pp 

1 1 1 


t S 2 


(4.57) 


(1  -t-  p)w^  = p,Pp^ 


t S 2 


(4.58) 


The  solution  ol  (4.57)  with  initial  conditions  w^(0)  = 0 and  w^(0)  _ 0 
is 

t 


The  solution  ol  (1.58)  is 


w = w (2)  . (t  - 2)w  (2)  . 

11  1 1+U 


(T))di]  (It 


1 2 2 (1.60) 


The  romainins  eciuations  ol  the  set  (1.53),  with  ^ given  by 

(1.55) ,  arc  solved  by  numerical  methods  '"the  translational  mode  n = 1 
may  also  be  solved  numerically  instead  ol  by  (1.59)  and  (1.60)'. 

The  dimensionless  radiallv  inward  displacement  w of  (1.55)  is  the 

o 

same  as  the  hoop  compressive  strain  .Vn  approximation  to  the  time  at 

which  the  elastic  del'ormntion  phase  ends  is  that  value  t at  which  w of 

o 

(1.55)  equals  the  compressive  vield  strain  e ol  the  shell  material; 

y 

expressefl  ma  t hema  t ica  lly 


w ( t ) = w 

o o 


(1.61) 


.\t  time  t,  the  radiallv  inward  velocity,  w (t)  = v , in  the  hoop  mode  is, 

o o 

bv  differentiation  ol  (1.55) 


I 


( 1 -t)  So( t --) 1 o 

- s,^e  “ |p^^(')d’ 


(1  .62) 


with  s^  = (o  - u)/2,  s^  = (v  ■ vi)/2,  and  v = s^  - s^  = - i(k^  - 2)'“ 

For  the  buckling  modes,  equations  (1.53)  are  solved  numerically  because  ol 

the  time  dependent  coellicienl  ol  w (an  analvtical  solution  is,  of  course' 

n 

po.ssible  it  w can  lie  neglected  in  t lie  coefficient), 
o 

fhe  initial  conditions  for  the  plastic  pliase  of  dotormation  arc 


w ( t ) = w 
n n 


w ( I ) = v 
n n 


n 2 1 (1.63) 


1.8  GOVER.MNG  EQL'.ATION  OF  THE  PL.\STIC  P1L\SE 

.-\s  n simplilied  description  of  the  complicated  plastic  response, 

buckling  is  taken  to  occur  in  accordance  with  the  Shanlcy  or  tangent 
(3 

modulus  theory.  In  this  theory,  it  is  supposed  that  the  memljrnno  and 
bonding  strains  combine  to  produce  no  strain  rate  reversal  on  the  shell 
cross  section.  Thus,  while  the  compressive  hoop  strain  increases,  the 
accompanying  cvirvaturc  increase  should  not  bo  excessive.  Experiments 
imply  conformity  with  this  condition  during  the  earliei'  [jart  of  the 
plastic  phase  because  low  buckling  modes  are  observed  (see  Figure  3.11). 
Strain  rate  reversal  does  occur  in  the  later  part  of  the  plastic  pliaso 
because  eventually  large  curvatures  develop.  Use  of  the  tangent  modulus 
theory  leads  to  an  overpredict ion  of  the  buckling  amplitude  because  the 
resistive  moment  is  much  higher  when  strain  rate  reversal  occurs;  how- 
ever, this  is  not  considered  a serious  handicap  of  the  theory.  Tiie 
advantage  of  tlie  theory  is  that  it  allows  a simple  analytical  treatment 
loading  to  conservative  results  from  a defensive  viewpoint. 

The  theoretical  development  lici'c  is  also  based  on  geometrical 

linearity,  that  is,  relatively  small  displacements  and  strain  (<  10 

percent).  To  obtain  the  simplest  possittlc  tlieory,  we  use  a constant 

tangent  modulus  E ; a variable  E can  reaflilv  ')c  incorporated  it 
h It 

wa  rrant  ed . 

The  stress-strain  curve  t Or  t tic  slicll  material  is  given  a liilinear 
representation  consisting  of  an  elastic  line  ot  slope  E and  a strain- 
hardening line  ol  slope  E , liaving  a point  ol  inti'rsect  inn  with  coordi- 

li 

nates  e and  , whicli  define  the  vield  sti-ain  and  vield  stress.  l?e- 
V V 

cause  materials  ol  interest  have  relativelv  low  values  ol  E , t lie  mem- 

li 

lirane  tlirust  is  considered  constant  and  e(|ual  to 

S k It-  k ilEe  = k HEw  a ( t . G 1 ) 

a \ a y a n 

in  terms  ol  plivsical  qu.intities. 

.">6 


According  to  the  tangent  modulus  tlicory  the  momen  t -cu  rvat  ui'o  change 
relationship  is 


M = E Ik 
It 


(1.65) 


By  lollowins  the  derivation  ol  the  K^vernint;  ec(uation  (d.l'i)  lor  the 
elastic  piiasc,  hut  usinn'  (1.65)  and  (1.61)  instead  o(  (1.1)  and  (1.8),  we 
find  the  dimensionless  lorm  ot  the  Kovernins  ecpiation  lor  the  iilastic 
phase  to  be 


^ . k k.k./^ 


■1  2 . 
w 9 


k w 1 


w|  = „p  (1.66) 


or 

where  t > t,  p _ Pp  + , and  k = E /K ; the  dimensionless  ((uantities 

3 h 

are  listed  in  Table  1.1.  The  initial  conditions  arc 


w(f),t)  = w(0) 


-(H  ,t  ) = v(H) 


(1.67) 


determinable  li'om  the  values  at  the  end  ot  the  elastic  phase. 


1.9  SYSTEM  OE  GOVERNlNCi  E()rATIONS  Ol  PUASTIt'  PHASE 

.Alter  substituting:  the  series  representations  (1.25),  (1.20),  and 
o r 

(1.2H)  lor  w,  p , and  p in  the  Koverninp;  equation  (1.12)  and  equating 


to  zero  the  coellicients  for  each  cosn  , we  obtain  the  lollowinu  svstem 
ot  equations; 


w + k w 
o 2 o 


= - r’, 


( 1 . 68  ) 


2 2 - or 

w 4 k (n  - l)(k  k n - w )w  = p,(Pp  - P ) t ^ I,  n ' 1 (1.69) 

n 2 13  on  n n 
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To  obtain  (1.68)  it  is  ajtain  assumed  that  w >-'<  1.  Introducinii  in  ('1.68) 

o 

and  (-1.69)  the  radiation  pressure  cootiicients  (l.-ll)  and  (1.-13)  results 
in 

w + .J.W  - (l,,'2)w  = - k w t > t (-1.70) 

o o o ()  2 o 


•J.V  )w  + ui3  w + 
n n n n 


2 2 - 
k (n"'  - l)(k  k n“  - w ) + uQ; 
2 1 3 o n 


w = uPp 
n n 


t S t n S 1 (‘1.71) 


To  obtain  (1.70),  wo  neglect  the  last  tenn  ol  (4. -11)  (this  term  is  retained 
in  the  summary  ol  Appendix  A).  The  accompanying  initial  conditions  are 

w (t)  = w w (F)  = V/  n > 0 (1.72) 

n n n n 

given  by  (4.61)  and  (1.62)  when  n = 0 and  (1.63)  wlten  n a 1. 


•1.10  SOLUTION  OF  GOVERNING  EQUATIONS  OF  PLASTIC  PHASE 

The  solution  ol  (1.70)  lor  the  hoop  mode  with  initial  condtions 
(1.72)  is 


w (-) 
o 


V f w ( a,  * 
o o 


k s ) 
2 1 


T 

1 

/V 


V • w ( g, 
o o 


s 


- e 


•’o<’ 


T))dT) 


(1.73) 


- t-, 


c 


where  t 


t 


W-)  ' 


-(  + .*)  2, 


f 


The  tifst  equation  ot  the  set  (1.71)  j;ovcrnint>  shell  translation 

(1  . ♦ “^l“l  ~ *^^’’^1  (1.7-1) 


and  lT)rmulas  (l.l-l)  t;ive  the  eoelticients 


q:  = - - 
1 2 


^ ^ 


V = 0 
1 


lor  t <.  2 


V = 1 tor  t 

1 


a = 0 B = 0 V = 1 tor  t 2 

1 1 

Hence  equation  (1.74)  is  replaced  by  the  tollowinp;  two  equations 


*1  ^ ‘*'^1  ~ 


(1-1  „>w  = _Pp 

1 1 


(-1.75) 


(-1.76) 


Equations  (-1.75)  and  (4,76)  arc  the  same  as  (4.57)  and  (-1.58)  so  that 
tlic  solutions  (4.59)  and  (-1.60)  extend  into  the  plastic  phase.  Physi- 
cally, this  result  simply  expresses  that  ritiid-body  translati<in  is  beinp; 
described  and  consequently  is  independent  ot  tlie  elastic  and  plastic 
de  tormat ions . 

Each  equation  ot  the  set  (4.71)  has  coefticients  that  arc  constant 
witliin  time  intervals  accordint;  to  (4.43)  and  (-1.44).  Thus 


n 

n 

- 

a - 0,  0 = 

n 

0,  V = 1/n  for 
n 

tat, 

A) 

and  scvc'ral  cases  can 

iiri  so  . 
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r 


For  t ^ t 2 'n,  (1/71)  becomes 


W + LiW  + 

n n 


k (n  - l)(k  k,  - w ) - an,  2 
2 1 3 o 


w = ,,Pp 

n n 


n S 2 (1.77) 


which,  for  initial  conditions  (d.72),  has  the  solution 


w (t)  = 

V + w ( s + a )1 

S T 1 

c /v  - 1 

V w (s  + a) 

n 1 

n n 1 ^ 

[_  n n 2 

o /v 


s^(t-?7)  s^(T-rj) 


p (t  i T])dT) 
n 


(1.78) 


where  t = t - t,  = (a  - v),'2,  s^  = -(p,  + v)  ' 2 , v = s^  - = 

, 1/2 


2 2 - 2 2 
+ l[an/2  + k_^(n  - l)w  - 


For  t S t and  t S 2/n  , ('1.71)  becomes 


..2  2 - o 

(1  ' a/n)w  + k (n  - l)(k  k,  n - w )w  = aPp  n 5 2 (1.79) 


1 3 


o n n 


which,  lor  initial  conditions  (1.72),  has  the  solution  (1.78)  with 

2 - 2 12 
s=-s  =^k(n  -l)(w  -kkn)(l+un)^  ,v=s  -s.  II 

122  ol3  " 12 

t < 2.  n the  initial  conditions  ot  (1.79)  are  w and  w at  t = 2/n  ob- 

n n 

tainable  from  (1.78)  with  t = (2  n)  - t. 


1.11  I INAL  DKFOKMXT ION 

The  final  dolormation  is  taken  as  the  shape  when  the  value  ol  the 

hoop  mode  dellection  determined  by  (1.73)  reaches  a maximum,  that  is, 

when  w (t  ) = 0.  The  mode  amplitudes  w are  then  determined  from  (1.78), 
o 1 n 

and  the  modes  are  superposed  to  produce  the  I inal  shape. 


(it) 


Table  T.l  (concluded) 


L Stiffener  spacing  (Figure  4.2) 

i Stiffener  thickness  (Figure  4.2) 

t Time 

0 Fluid  density 

0 Shell  material  density 

s 

1/2 

c = (k/p)  Fluid  sound  speed 

1 2 

c = (E/p  ) Shell  sound  speed 

s s 

K Fluid  bulk  modulus 

E Shell  modulus  of  elasticity 

E Shell  hardening  modulus 

h 

p Fluid  pressure  on  shell 

k,  = + X - 1)  - (>  - i)(:  - i)‘]/\:^(:  + \ - d 

k = (c  V X - i)/x: 
a 


1 
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TI I KOKK T I C A L KS U LTS 


5.1  NLAIERICAL  DATA 

An  o.xternally  slillonecl  model  cylindrical  shell  (Figure  1.2a)  has 
been  chosen  to  illustrate  the  results  predicted  by  the  theory  described 
in  Section  1.  The  data  lor  this  numei’ical  example  are  listed  in  Table 

5.1. 

The  model  was  chosen  to  turnish  realistic  values  ol  tlie  reciuirod 

input  constants  k^,  k^ , 4,  and  In  this  example  the  value  ot 

E was  chosen  as  an  average  value  lor  steels  over  plastic  strains  up 
h 

to  about  live  times  the  yield  strain;  it  was  fdund  that  linal  delorma- 

tiftns  were  insensitive  to  E lor  k in  the  range  0.02  k <0.1.  With 

h 5 d 

geometric  scaling  and  the  same  material,  the  values  ol  tlu'se  constants 
do  not  change,  so  the  results  are  applicable  to  a lull-scale  slillened 
shell.  Similarly,  the  values  ol  P and  T arc  realistic  and  independent 
o I scale  . 

5.2  XU.\IEI!  ICAI,  IlESFETS 

.\  solution  was  obtained  in  accordance'  with  t lie  treatment  outlined 
in  .Appendix  B.  Various  illustrative  portions  ol  the  solulion  are  shown 
in  Figures  5.1  through  5.10.  For  this  example  it  was  lound  that  Iti  motle 
provided  an  adequate  desci’iption  ol  the  shell  ri-sponse  becnusi-  higher 
modes  are  not  amplitied  during  the  motion  (trigonometric  oscillation 
instead  ol  hyperbolic  growth).  The  results  are  discussed  in  terms  ol 
d imension less  quantities. 

Figures  5.1  to  5.(1  show  t lie  growth  ol  the  displacement  and  velocitv 
amplitudes  lor  morlcs  n O to  n = 5.  Figiii'c  5.1  shows  the  tuhavior  ol 

(w 


t lio  hoop  mode.  The  associated  luuncrical  results  show  that  plastic  1 low 
occurs  at  t = I = 1 . 15  wi  t li  w = w = 0.002  and  that  tlie  lioop  motion 
ceases  at  t ^ = 1.13  (incident  pulse  duration  is  T = 3),  at  which  time 
w = t).00«7.  Tile  maxiiiuun  hoop  strain  is  about  -1.2  e . The  plot  lor 

t)  y 

the  velocity  shows  an  initial  rapid  increase  lollowed  by  a gradual  de- 
crea.se  itntil  plastic  t low  occui’s,  whereupon  the  velocity  increases  until 
pressure  tmloadint;  commences  at  T = 3. 

1- ieur<'  5.2  shows  tlu’  liehavior  ol  the  translation  mode.  The  maxi- 
iUm  displacement  . e-  a riuiil  lioilv  is  w^  - O.OIH,  which  is  1 .8  percent 
ol  the  radiLii  o'  l.i-  wall  :i  i km  ■- -si 's  . The  translational  velocity 
o\orshoot  ■ the  :■  .'ll  pill  ■ ■ I'ti'li-  I'  locitv  ol  0.005  to  reach 

'w  ^ o.iiii()i.  • ' • . . M 1 tends  to  the  incident  wave  parti- 
cle i . 1 o(  I I 'iin,  oi-ciirs  at  T = 3.  Some  reverse 

' rai’  ■ 1 i i . I < j-  i _ 30.  Computations  carried 

oil'  ''T  ■ ■ sh-  w that  the  shell  comes  to  rest. 

■' 1 s ' , i ' ■ ■ I'o  ' >r  1 ii  .or  times  with  no  unloadint;  show 

Iha'  til  ;i  ; I : ■ il  a i t v I'lentuallv  assumes  the  incident 

t 0|)  W O '’ll'  . 1 ' . 

1-  1 . u 1 1 • -w  - 'III  di  • I o I"  a t 1 on  amplt  tilde  in  the  second  mode, 

which  IS  ’hi  I .iliie,  o I l•.2d  loiirier  component  . The  maximum  dis[)lace- 
nii'iit  ol  w o.iMil’  occiii's  at  I mil  t r 2.55  and  is  about  one-hall  ol 

V 

I III  ( . ill  'loop  . I(li  d i P 1 .11  1 m t o I 0.0087. 

The  behavior  ol  the  i rails  1 a t i ona  1 motion  ilepends  on  the  unit  resistance 
liinction  (t  ) chosep.^to  approximate  the  exact  lunction.  In  this 
example  t (1)  e was  used  to  provide  a iinal  shell  translational 

velocity  ei^ia 1 to  the  step  incident  wave  particle  velocity.  11  the 
two  slraiijht  line  .ipprox  1 ma t i on  (1.39)  is  used,  the  final  shell 
translational  velocity  is  about  twice  the  (larticle  velocity.  I he 
mathematically  correct  velocity  probably  lies  closer  to  the  latter 
value.  This  i nvest  ip.it  i on  was  not  pursued  because  translational 
behavior  has  no  inlliieiu'e  over  the  Iinal  del  Drmat  i on . 


Figures  5.‘1,  5.5,  and  5.6  show  the  third,  lourth,  and  liith  mode 
displacement  amplitudes  witli  maxima  at  tlic  slicll  unloading  time  ol 
t = 1.13  dcteniiinod  by  the  hoop  mode.  These  maxima  arc  0.0025,  0.0021, 
and  0.0015,  respectively.  Corresponding  ligurcs  loi'  modes  6 through  10 
are  not  included  here,  but  they  sliow  continued  trends  of  d imin  i siting 
maximum  amplitudes  and  they  all  occur  at  or  near  t^. 

Figure  5.7  is  a plot  of  the  maximum  displacement  amplitudes  of 
modes  0 to  10  at  times  t,  T,  and  t^.  The  jtlot  allows  a comparison  ot 
modal  amplitude  increase.  For  modes  above  about  n = 7 the  gretwth  is 
negligible;  in  tact,  those  amitlittides  arc  predicted  by  the  theory  to 
be  oscillatory.  From  the  distribution  ol  figure  5.7  it  is  seen  that 
the  deforming  shape  of  the  shell  cross  section  cotttains  mostly  the 
Fourier  components  n = 0 ami  n = 2-6  (n  = 1 is  rigid  body  translation). 

Figure  5.8  shows  Cartesian  plots  of  the  deformed  shape  of  the 
shell  at  times  t,  T,  and  t^.  In  this  example,  the  shape  at  t (when 
plastic  flow  starts  in  the  hoop  mode)  is  generally  amplilied  during 
plastic  buckling.  However,  realistic  examples  exist  where  the  plastic 
buckling  has  substantial  higher  modes,  especially  for  lower  values  ol 
the  strain-hardening  modulus  (e.g.,  k^  < 0.02). 

Figure  5.9  shows  a polar  plot  of  the  final  deformed  shape.  The 

o 

maximum  radially  inward  displacement  occurs  at  H = i -It)  and  at  t tic 
o 

Itack  at  0 = 180  (see  also  Figure  5.8);  tlte  displacement  magnitude  at 
o 

y = i 40  is  0.0111,  wliich  is  1.66  times  tlie  lioop  mode  deflection 
w^^(t^).  Tliis  collapse  shape  lias  been  oltservcd  in  scale  model  tests 
carried  out  at  SKI  (see  Figure  3.11  lor  a qualitative  comparison). 

Figure  5.10  shows  the  outer  I liter  strain  around  the  shell  at  times 
t,  f , and  t^.  .At  time  t^  tlie  neutral  axis  hooit  strain  is  -0.0087,  and 
t lie  departures  front  this  value  are  t 0.0033.  Although  an  analysis  to 
examine  the  extent  of  strain-rate  reversal,  if  any,  has  not  vet  been 
included  in  the  theoretical  treatment,  the  strain  curves  imply  that 
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reversal  was  absent  in  tliis  example  (curves  do  not  cross  eacli  otlier). 

The  lower  the  hardening  modulus  value  , the  greater  the  possibility 
ol  strain-rate  reversal. 

Figures  5.11  through  5.14  show  modes  n = t)  to  3 lor  the  dil traction 
o r 

pressure  Pp  - p . Figure  5.11  shows  the  pressures  associated  witli  tlie 
n n 

hoop  motion.  The  ditlraction  pressure  monotonically  tends  to  the  Irec- 
lielil  pressure  of  P = 0.005  until  pressure  unloading  occurs  at  T = 3. 

The  radiation  pressure,  plotted  as  a pressure  relief  caused  by  the 
shell  motion,  is  seen  by  comparison  with  Figure  5.1  to  depend  primarily 
on  the  shell  velocity  during  the  initial  elastic  phase  of  deformation. 
During  the  plastic  phase,  the  pressure  relief  increases  in  a manner  simi- 
lar to  the  shell  velocity  but  at  a slower  rate  because  of  tlic  compensating 
effect  of  the  convergent  fluid  afterflow.  Dui'ing  the  clastic  phase,  tlie 
total  pressure  rises  to  a plateau  of  0.003,  whicli  is  maintained  until 
incident  pressure  unloading  at  T = 3.  After  time  T = 3,  tlie  total  pres- 
sure regains  a comparable  value  because  of  the  fluid  afterilow  and 
inertia . 

Figure  5.12  shows  the  pressures  associated  with  the  rigid-body 
translation  of  file  shell.  The  main  feature  is  tlie  relatively  small 
total  pressure  brought  about  by  the  pressure  relief  caused  by  the  shell 
mo  t i on . 

Figures  5.13  and  5.11  show  the  pressures  associated  with  shell  de- 
fomation  in  the  second  and  third  modes.  As  in  the  case  of  the  transla- 
tional first  mode,  the  pressure  rclie!  is  considerable.  Pressure  magni- 
fudes  generally  decrease  and  oscillation  Irequcncies  increase  as  the 

mode  number  increases.  i'hc  modal  amplitude  velocities  w corrospomi 

n 

primarily  to  the  diffiaction  pressures  in  the  elastic  phase. 

Figure  5.15  shows  the  dil  tract  ion  pi-cssiire  at  the  tront  ((1  = 01, 
top  and  bottom  (0  - - 2) , and  rear  (9  - ")  ol  a rigid,  I i xcd  cylinder 
siitjjected  to  the  incident  rectangular  pulse  ol  pressure  P --  0.005  and 


T 


duration  T = 3,  as  described  by  the  haniionics  n = 0 to  10.  On  arrival 
oi  the  wave  at  the  I ront  oI  the  cylinder,  the  incident  pressure  doubles 
in  amplitude  to  0.01,  which  is  in  accordance  with  acoustic  reflection 
from  a rigid  plane  wall.  Thereafter,  the  pressure  decays  to  the  incident 
pressure  because  of  dil tract  ion  around  the  cylindrical  surface.  .At  the 
top  and  bottom  of  the  cylindoi'  the  pressure  history  is  not  substantially 
different  from  the  incident  wave  history.  .At  t!ic  rear  of  tlie  cylinder 
the  pressure  rises  gradually  to  tlio  incident  pressure.  Unloading  demands 
the  sustaining  of  a negative  pressure  in  the  fluid,  which  physically 
moans  cavitation  witli  zero  pressure.  In  actual  underwater  explosions, 
tliis  effect  would  t)o  reduced  because  tlic  pulse  cutoff  would  not  be 
instantaneous . 

Figure  5.16  shows  the  total  pressure  at  the  front,  top  and  bottom, 
and  roar  of  a cylindrical  shell.  .At  the  front  the  pressure  rises  to  a 
value  of  0.009  (U  = 0.005)  tnit  is  rapidly  relieved  by  the  shell  wall 
motion;  the  pressure  at  0=0  drops  to  a value  less  tlian  0.001.  .Aftcj' 
the  minimum  pressure  is  reached,  the  pressures  at  0 = 0,  ± ”2,  and  ~ 
rise  smootfily  togetficr  to  a plateau  of  about  0.003.  .A  shoj-t  duration  of 
negative  pressure  is  predicted  immediately  following  unloading;  tliis 
condition  is  not  so  severe  as  tliat  shown  in  Figure  5.15  lor  the  rigid, 
fixed  cylinder.  It  is  probatile  that  a realistic  pulse  cutoff  will  re- 
move tills  unrealistic  feature. 

.A  striking  qualitative  comparison  witli  experiment  is  afforded  by 
Figures  A-3  and  3.12.  Theory  and  experiment  fiot h show  a sharp  rise 
followed  by  a sharp  decay  and  gradual  rise.  The  final  decay  differs 
because  the  theory  is  based  on  an  idealized  rectangular  pulse.  .Aside 
from  the  final  decay,  the  main  difference  in  the  pulses  is  the  pressure 
level  in  the  valley  between  the  sharp  decay  and  gradual  rise.  This 
pressure  level  depends  on  the  nature  of  the  structure,  ami  since  the 
structures  dilfer,  these  levels  sliould  not  be  t lie  same. 


I'itiurcs  5.17  to  5.21  exhibit  the  liislory  o 1 the  total  pressure 
distribution  around  the  cylindrical  shell  by  moans  of  plots  at  times 
t = 0.5,  1.15,  2.0,  3.0,  and  1.13.  uniform  pressure  distribution  is 

almost  established  at  time  t = t = 1 . 15  (I'iKurc  5.18)  oven  t liough  at 
this  time  the  inciilent  wave  front  has  travelled  only  1.15  shell  radii. 

•\  t the  shell  engulfment  time  t = 2,  the  pressure  distribution  is  unilonii 
(Fixture  5.19)  and  remains  uniform  at  least  until  the  pulse  duration  time 
t = T = 3 (Figure  5.20).  Departure  from  unifonnity  is  noticeable  at 
time  t = t^  = 4.13  (Figure  5.21).  These  results  arc  consistent  with 
the  total  pressure  results  displayed  in  Figures  5.11  to  5.14,  whicli  show 
that  the  hoop  mode  Fourier  component  ot  pressure  rapidly  roaches  a con- 
stant value  and  is  much  greater  than  the  remaining  components. 


Table  5.1 


DATA  FOR  STIFFENKD  CVI.INDRICAI.  SItKI-L 


Proper  I y 


Symbol 


\aliio.s 


Internal  radius 


1.11  in 


3 . 66  cm 


Wall  thickness 


0.0167  in 


0.0-1 21  cm 


Still  ener  dep  1 li 


0.09H5  in 


0 . 25  cm 


Stilt  cue  r width 


0.030  in 


0.0762  cm 


Stillener  spacinn 


0.135  in 


1 . 105  cm 


Shell  density  (steel) 


7.8  si/cm 


Fluid  density  (water) 


1.0  i;  cm 


Shell  clastic  modulus 


30  X to  psi  206,897  MPa 


1 1 u i d Ini  1 k modu  1 us 


0.33  X 10  psi  2,276  MPa 


Yield  St  ress 


60,000  psi 


Yield  St  ra i n 


3 = C'  E 

y y 


Shell  liardenitu;  modulus 


3 X 10  ()Si 


20,690  MPa 


Stillener  wiilth  parameter  X = L/i 


Stillener  di-ptli  parameter  = II  h 


\ rea  Coe  t I i c i en 1 


0 . 221 5 


Second  area  monii-n  t 
l lie  I I i c lent 


0. 167  1 


Yi'Utr.il  axis  coellicient 


0.11  15 


Nell  I ra  1 a X I 'll  ili  p t h 


/ k II  2 

o / 


0.020  I 


i Shi ' 1 I r.id  I us 


a a . 7. 


Values 


1,509  m/scc 
5.151  m/sec 

11.655 
0.1 
8.  165 


0.002 

0.005 

1650  psi 

11.38 

MPa 

925  psi 

6.38 

MPa 

3 


75  u.sec 


FIGURE  5.4  DISPLACEMENT  AND  VELOCITY  AMPLITUDES  FOR  THE 
THIRD  MODE  In  3) 
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FIGURE  5.13  DIFFRACTION.  RADIATION.  AND  TOTAL  PRESSURE 
ASSOCIATED  WITH  SECOND  MODE  (n  21 


FIGURE  5 14 


DIFFRACTION.  RADIATION.  AND  lOTAL  PRESSURES 
ASSOCIATED  WITH  THIRD  MODE  (n  3) 


DIFFRACTION  PRESSURt  AT  FRONT,  TOP  AND  BOTTOM 
AND  REAR  OF  RIGID,  FIXED  CYLINDER 


FIGURE  5 15 


TOTAL  PRESSURE  AT  FRONT,  TOP  AND  BOTTOM,  AND 
REAR  OF  CYLINDRICAL  SHELL 


FIGURE  516 
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FIGURE  5.17  TOTAL  PRESSURE  DISTRIBUTION  AT 
TIME  t 0.5 
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FIGURE  5.18  TOTAL  PRESSURE  DISTRIBUTION  AT  TIME 


6.  Pl'ESSTOE- IMPULSE  ISODAMAGE  KEPUESENTATIOX 
AN'D  COMPARISON  WITH  EXCESS  IMPULSE  RULE 


An  ollcctive  way  oC  representing  the  damage  to  a structure  by  pulse 
loading  is  by  an  isodamage  curve  in  the  pressure-impulse  plane.*  For 
an  incident  rectangular  pulse  of  pressure  P and  duration  T,  the  impulse 
is  I = PT . Damage  to  a submarine  shell  structure  may  be  taken  as  t lie 
maximum  pemanent  radial  displacement  6;  in  dimensionless  tenns  this 
displacement  may  be  represented  as  a traction  ot  the  shell  radius  or  a 
multiple  ot  the  shell  thickness.  A pressure-impulse  curve  is  therefore 
the  curve  in  the  P-I  plane  on  wliich  5 is  a constant;  the  curve  connects 
rectangular  pulses  ot  equal  c t tcct i vcncss  that  range  from  ideal  impulses 
to  step  loading. 


6.1  PI  CURVES  FOR  A CYLINDRICAL  SIKLL 


Tltc  thcorv  dcscrilted  in  Section  I was  used  to  generate  P-I  curves 
for  an  externally  st  i I tened  cylindrical  shell  having  t tie  properties 
listed  in  Table  5.1.  .Additional  cl ar i 1 ica t ion  of  dimensional  del  ini t ions 


The  pressure  and  duration  ot  t lie  incident  rectangular  pulses  are 
represented  by  P and  T and  the  damage  is  represented  nond imen si ona 1 1 y by 


' = maxlw(u,t  ) - w (t  ) cost) I 

a ' I It 


(6.1) 


that  is,  by  the  maximum  inward  radial  displacement  at  the  t ime  when 

the  hoop  mode  inward  displacement  w^^(t)  reaches  a maximum;  t^  is  given 

bv  w (t  ) 1:0  and  the  rigid  bodv  translat  ional  mode  is  siditracted  I rom 
I)  t 


«6 


the  total.  Uimcnsional  Iv , is  the  maximum  linal  inward  di  splacomont 

a 

divided  by  the  shell  radius  a. 

Pressure-time  (P-T)  eui'ves  tor  isodama^e  v'alues  ot  ^ equal  lo  0.02 

a 

0.03,  and  0.04  are  shown  in  I'ifiure  6.1  and  the  corrospondin}>'  P-1  curves 
are  shown  in  Figure  6.2.  The  overall  shapes  ol  llio  curves  resemble  rec- 
tangular hyperl)olas  and  arc  typical  (5i  critical  load  curves  lor  many 
structures.  The  horizontal  asymptote  P = 0.0025  gives  the  pressure  in 
a step  incident  wave  that  is  required  to  .jvtst  cause  plastic  dclormation. 
Cross  plots  that  show  how  damage  gradients  steepen  with  increasing  T and 
P arc  shown  in  Figures  6.3  and  6.4  (compare  Figures  3.13  and  6.4). 

6.2  PIIYSIC.AL  QUAN'TITIFS 

The  f values  ot  the  isodamage  curves  ot  Figures  6.1  and  6.2  arc 
a 

the  ratios  ol  ma.ximum  permanent  inward  displacement  to  t lie  original 

radius;  ttiai  is,  the  damage  displacement  curves  represent  2,  3,  and  4 

percefit  ot  t lie  original  radius.  Because  tlie  ratio  of  radius  to  wall 

thickness  is  a h = HH,  tlie  isf)damagc  curves,  as  multiples  ot  wall  tliick- 

ness,  are  = 1.8,  2.6,  and  3.5.  I'n  i t values  ot  P = 0,001,  T = 1,  and 
ii 

1 = 0.01  correspond  to  an  incident  |)ressure  ot  p.  = KP  = 330  psi  , a 

radius  transit  time  ot  t = aT/c  = 25  ,;scc,  and  an  impulse  ol  82.5  psi- 

d 

msec.  For  an  acoustic  stc>p  wave  traveling  lioii  zont  al  ly  in  water,  the 

I rec  surlace  cuts  oil  the  pulse  experienced  al  a depth  11  at  a time  equal 

to  tlie  transit  time  t = 11/ c. 

d 

For  t lU'  model  cylindrical  sliell  stdijected  to  a rectangular  incident 

pulse  with  a duration  ef|ual  to  a transit  time  ol  two  diameters  (t  = 

d 

100  sec,  T _ 4),  a pressui'e  o|  p = 2046  psi  (P  = 0.0062)  will  cause  a 

i 

permanent  dclormation  ol  3.5  wall  thicknesses  (*  = 0.04);  the  corre- 

a 

sponding  incident  impulse  is  204  [isi-msec  (I  = 0.0248),  The  static 
vield  pressure  is  p = ^ (k  H)  a = 925  psi,  giving  a pressure  I'atio 


kL 


p./p  = 2.21;  lluis  the  excess  impulse  is  I = (p  - p )t  - 112  psi- 
1 y X i y d 

msec  (1  = 0 . 0136) . 

As  another  example,  il  the  pulse  has  a transit  time  oi  lour  diameters 

(t  = 200  usee,  T = 8) , a pressure  ol  p = 1250  psi  (P  = 0.0083)  will 
d i 

also  cause  a pemanent  delormation  ol  3.5  wall  thicknesses  = 0.01); 

a 

the  pressure  ratio  is  now  p,  p = 1.35.  The  incident  impulse  is  250 

1 y 

psi -msec  (I  = 0.301)  and  the  excess  impulse  is  I =65  psi -msec 

X 

(I  = 0.0079). 


6.3  COMPARISON  WITH  EXCESS  IMPULSE  RULE 

■\n  examination  ol  the  isodamape  curves  shows  that  the  excess  im- 
pulse depends  on  the  P-I  or  P-T  point  selected.  For  example,  along 

curve  ' = 0.01  ol  Figure  6.1,  the  cl  intension  less  excess  impulse  rasiges 

a 

Irom  about  O.OIT  to  0.008,  a di 1 lerence  ol  almost  a lactor  ol  two,  as 
P-T  points  range  Irom  (0.007,  3.1)  to  (0.0038,  8.0).  Hence,  excess 
impulse  is  not  constant  lOr  cc^istant  damage,  and  therelore  excess  im- 
pulse is  not  as  accurate  a measure  ol  tlic  damage  potential  ol  undci-watei- 
shock  waves  as  the  P-I  representat ion . 
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FIGURE  6.3  DAMAGE-DUhATION  GRADIENT  CURVES 
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Appendix  A 


KXC’KSS  IMPULSK  RULE 

liu'ure  A.l  shows  a skctcli  ol  a simplilicd  li’co-l'ield  pressuro-t  imo 

curve  1 rom  an  underwater  explosion.  Here  i>  is  the  static  in-essurc  that 

c 

will  collapse  tlie  hull,  and  p^  is  the  local  hydrostatic  pressure.  The 
excess  impulse  is  the  shaded  area.  The  Excess  Impulse  Rule  states  tiiat 
hull  damaf;c  depends  on  the  magnitude  ol  t lie  excess  impulse. 

Parnell  and  Schauer  introduce  i lie  Excess  Impulse  Rule  as  follows. 
The  submarine  hull  is  regarded  as  a risid-plast ic  system  ol  one  dettree 
ol  freedom.  .As  indicated  in  lip;ui’c  A. 2(a),  the  loati  on  the  hull  is 
taken  to  be  uniform,  and  the  hull  is  assumed  to  deform  uniformly  inward. 
The  hull  material  is  taken  to  liave  the  rin'id-plast  ic  st  ress-s  t I’a  in  curve 
siiown  in  1'ip.urc  A. 2(b).  Tills  representation  ol  a submarine  hull  can  be 
expected  to  be  valid  foi’  distributed  loads  for  which  the  delormations 
greatly  exceed  the  elastic  rant;e. 

The  equation  of  motion  ol  the  model  ol  Fi;;ure  A. 2 is 

mv  = i)(l)  - P (A.l) 

o 

where  m is  the  mass  per  imit  ai'ea  of  the  hull,  y is  the  inwai'd  acceici'a- 

tion,  p(  I ) is  the  pressure  act  inp;  on  the  hull,  and  P is  the  pressure 

o 


* 


U.  (■ . Parnell  and  11.  M.  Schauer,  "Submarine  Hull  Rc-sponse  and  Damage 
Development,"  Chapter  (5  in  Handbook  of  Nuclear  Explosions:  Part  II 

Ellects,  Report  DASA-1210  II  ((f)  andC-21HS  (.lunc  IlUit;)  . 


tliat  must  be  exceeded  to  initiate  deformation.  Denotinji,  as  beloi’e, 

the  static  pressure  reciuii'ed  to  collapse  the  hull  by  p and  the  local 

hydrostatic  prcssui'e  by  [>  , we  haye 

h 


I’ 


o 


p 

C 


P 


h 


and  (A.l)  becomes 


my  = p(  t ) - ( p ■ P,  ) ( A.  2) 

c h 

To  solyo  (A. 2)  in  terms  ol  t lie  Iree-ticld  pressure,  we  must  know 
the  relationship  between  the  tree-1 ield  pressure  and  p(t),  the  load 
act  int;  on  the  hull.  for  this  relationship,  Parnell  ami  Schauer  uive  , 
at  the  point  ot  normal  incidence  during  the  enpullment  iihase 

p(  t ) = 2 p(  t ) - n cy  ( A . 2 ) 

whore  p(t)  is  the  local  Iree-tield  pressure,  p and  c are  the  density  and 
sound  speed  ol  water,  and  y is  the  inwai'd  velocity.  Thus,  the  trec-licld 
pressure  is  doutiled  by  rcl  lection  and  reduced  by  the  term  pey  due  to  hull 
motion.  The  pressure  calculated  usinp;  this  interact  ion  law  is  compared 
with  the  measured  Irce-field  pressure  and  the  measured  pressure  on  the 
surlace  ol  a model  toi'  a typical  example  in  fipure  A. 2.  The  calculated 
[iressure  history  reproduces  the  initial  spike  reasonably  well , but  repre- 
sents the  remaituler  rather  poorly. 

for  the  pos  t enttti  1 I men  I phase,  Parnell  and  Schauer  take',  inste'ad  ot 
(A. 2)  , 


p(  t ) 


p(  t ) 


pey 


( A . n 


lor  the  intoraclion  law.  UsinK  (A.l)  in  (.\.2).  Ihcy  liiui  as  an  approxi- 
mate relation  between  the  inward  displacement  d and  the  excess  impulse 

1 , 

X 


ccd  = 1 (A. 5) 

X 

Because  ot  the  many  approximations  involved  in  deriving!  (A. 5),  Painell 
and  Schauer  abondon  the  theoretical  approacli  at  this  point  and  a(\opt  the 
functional  loafl-dc  fleet  i on  relation 


d = 1(1  ,11) 

X 


(A.B) 


where  f is  a function  to  lie  determined  from  ex  '>  -ital  data  and  H repre- 
sents a combination  of  parameters  charactoriz  ' hull  design  and  depth 
of  submergence. 

To  determine  (A. 6),  Pamcll  and  Schauer  introduce  the  two  parameters 


where 

d = maximum  stil  loner  dot  loot  ion  (indies) 
li  hull  plat  inf;  thickness  (indies) 

,S  still  one  r 1 ac  t o i' 

s rclerence  stilliMier  factor  (S  U.olB) 

o «' 

vii'ld  stress  ot  liull  material  (psi) 


n lennce  yn  ld  .stress  (j 


, SI), (ton  |isi  ) 


and  the  stii'fener  factor  is 


This  is  the  K/coss  Impulse  Uule.  Thus,  the  Excess  Impulse  Uule  is  an 
empirical  relation  based  on  the  correlation  ol  lest  results  indicated 
in  I'iRurc  A.l. 
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TIME  t 


MA  in  1 5 


FIGURE  A 1 DEFINITION  OF  EXCESS  IMPULSE 
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Appendix  B 


SU.MMAriV  OF  .UATIiL’.MATICAL  TREAT.\[F.\T 

This  appendix  collects  the  goveiTiinK  equations  and  solutions  lor 
dynamic  overall  buckling  by  a rectangular  incident  pulse.  .Also  in- 
cluded are  some  improvements  to  the  approximate  results  given  in  Section 
4 that  were  suggested  by  tlie  study  of  numerical  examples.  The  treatment 
given  ■'  this  appendix  was  used  to  generate  the  numerical  results  of 
Sc'  The  description  is  developed  by  modes. 

. Linuamental  Mode  (n  = 0) 

The  governing  equations  for  elastic  (0  < t < t ) , and  plastic 

(t  > t)  motion  in  the  hoop  mode  are 

w + k w = u(Bp”  - p^)  0 < t < t (R.l) 

o 2 o 0 0 

w = u(Pp  -p)-kw  t>t  (B.2) 

o (1  0 2 o 

Expressions  for  t lie  diflraction  and  radiation  pressures  are  (1.40) 

and  (4.41),  oljtained  by  substituting  the  approximate  unit  resistance 

tunction  (4.38)  in  lormulas  (4.3(5)  and  (1.33).  For  stof)  loading  a bolter 

* 

appi’oximat  ion  to  I lie  diffraction  pressure  is  given  by  Mnev  and  Pertsev 
as 

o -O.Uilt  -1.37(5t  -8.()t 

p^^(t)  = f^^(t)  = 1 - 0.282O  - 0.391e  - ().327e  (B.3) 


Ye.  Mnev  and  A.  K.  Pertsev,  ’’Hydroelasticity  of  Shells,  ” English 
translation  of  Fon.'ign  Technology  Division,  \Vr  i gh  t -Po  1 1 e rson  Air 
Force  Base,  Ohio,  Report  AD-731  (5  1(5  (24  August  1971). 
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Thus,  loi'  a i-octanKular  pulse  ol  diu-ation  T, 


p^^(l)  = ll(t)  - H(t  - T)1  (1  - T) 


(B.  1) 


Equations  (B.l)  and  (B.2)  may  now  be  wi’itton  in  tlio  form 


o o 


t 

‘ ^*^0  “ + (u/da)  I w^^(t)o  “""’dT  = pPi)”(t) 


I 

1 

0 


0 ■'  t < t 


(B.5) 


t 


Wq  ‘ “ (m./2)w^^  + (u/la)  \ w_(t)g  = a.Pp''(l)  - k„w 


t > 1 


(B.(3) 


with  the  dil fraction  pressure  given  by  (B.4)  and  (B.B).  Equations  (B.5) 
and  (B.6)  are  solved  by  numerical  methods. 

Translational  Mode  (n  = 1) 

The  equation  governing  translational  motion  is 


o 1' 

- V^) 


(B.7) 


Expressions  lor  the  dillraetion  aiul  radiation  (jressure  are  (1.12) 
and  (1.13),  obtained  by  substituting  the  approximate  unit  resistance 
lunction  (1.39)  in  formulas  (1.36)  and  (1.33).  lor  step  loading  a betler 
approximat  ir>n  to  the  ditlraction  pressure  is  given  by  Mnev  and  Perl  sev 
as 


p‘^’((  ) - 1^(1  ) 


2.97 


-0 . 6161 
l'  CoS  ( ( 


3) 


-1 . 7251 
e 


+ 


./  ->.5t 

0.6  le 


-16.nt 

I 


(B.B) 


ini 


L 


p^(t)  = II(t)  - li(t  - T)f^(t  - T) 


(B.9) 


Anain  for  step  loadiiifi',  a better  approximation  to  the  radiation  pressure 


p^  is  obtained  by  introducing  the  approximate  unit  resistance  function 
-t  2 

i/)^(t)  = e instead  of  (4.43)  in  formula  (4.33).  This  function  gives 


a final  shell  translational  velocity  equal  to  the  step  incident  wav'e  par' 
tide  velocity.  The  resulting  governing  equation  is 

t 


f -(t-T)/2  o 

- (u,  2)w^  ^ (p,/“l)  I w^(t)o  dr  = pPp^(t)  (B.iO) 


0 


with  the  diffraction  pressure  given  by  (B.8)  and  (B.9).  Equation  (B.IO) 
is  solv'ed  l)v  numerical  methods. 


Buckling  Modes  (n  > 2) 


The  governing  equations  for  elastic  (0  < t < t)  and  plastic  (t  > t) 
motion  in  the  higher  modes  are 


w . k ( n - 1 ) 
n 2 


2 1 o 1' 

k n - w (t)  w = „(I’p  - p ) O < t < t 


o j n 


(B.ll) 


2 2 - o r 

w +k(n  - l)(kk,  n - w)w  = .*(Pli  - P ) 


1 3 


t > t (B. 12) 


with  tile  diffraction  presstirc  p , obtained  bv  stibst  i t ut  ing  (1.39)  in 

n 


( 1 . 3 (i ) , is 


p ( t ) a ( t ) - b ( t ) - ( n 

n n n 


2)j"  b^^(')d'  - 


t-2  n 

(n  2)11(  t - 2 n)  I b (-)d-  (B.  13) 

J ” 

0 


in  which,  lor  a ttnil  rectangular  incident  pulse  ol  duration  T, 


lO.'i 


J 


Unloading  (t  > T)  (Cases  lor  T > 1): 

t 

n 

-2/n 


= I (I  ) 
n 


K (t)  . 
n 2 


1 (t  • 
n 


H U 
n 


n 

T)  - - 
o 


(-  - T)di 


max( ( -2 . n ,T) 


-I  K (T)dT 


n 
J 

t-2.  n 


1 (t  - T) 
n 


4 


K (t 


n 


t 


max( t-2  n ,T) 


2 . - < I T - 2,  p"  =-I(t-T)  + K(t-T) 
n n n n 


-tjv 


(t  - T)dT 
I 

niax(  t -2,'n  ,T) 


O 11  f 

’ = - o 


T . 2 

(T  - T)d" 


T . 2 < t < T 4 2 + 


Tlic  radiation  pressure,  siivon  Ity  (1.13)  and  (-I. 11),  is  substituted 


in  the  ftovci'iiiny  oquations  to  produce  the  followiiift  two  equations  itjr 


elastic  (()  t t)  and  plastic  (t  > t)  motion; 


...  r 2 2 -|  o 

(1  ^ uv  )w  + w j k (n  - l)(k  n - w ) + Ltd  w = ul’p  (B.15) 

n n n L 2 1 o 


. r 2 2 - 1 o 

(1  * U.V  )w  + e/3  W + k (n  - l)(k  k n - w ) • eO!  w = I’p  (B.K5) 

n n '■^n  n [_  2 13  o nj  n n 


with  the  dillraction  |)rcssure  p (>ivcn  by  (B.13)  and  (B.ll),  and  t lie  values 

n 

of  the  parameters  a , n , \ niven  by  (l.ll),  that  is 
n ^n  n 


a = -n/2  Q = 1 

n 


a = 0 
n 


V = 0 (or  t ■'  2/n 
n 


V = 1 n tor  t > 2/n 
n 


In  Koncral,  equation  (B.15)  KovorninM;  elastic  motion  is  solved  numer- 
ically because  ot  the  term  w (t).  Ilowovci’,  for  practical  values  of  the 

o 

parameters  k , k , k , and  e,  and  fOr  n sufficiently  lar/;e,  w (t)  mav 
12  3 o' 

be  ncalected.  For  Utesc  cases,  solutions  are 

I 


w 

n 


<>  , -|J.(t-T)  2 , 

p (t)c  sinv  (t 

n n 


T)dT 


(B.17) 


0 < t < 2 n 


where  v 

n 


(u./2) 


1 '2 


2 2 2 
\ = k k (n  - ])n  - nii/2 

n 12 


lOH 


w = w cos\  (t  - t)  + — sin'  (t  - t) 
n n n ' n 


(1  • li/n) 


t 

f " 

— \ p (T)sin' 

' J n n 


(t  - T)dT 


- < t < t (B. IS) 
n 


2 2 
k k ( n - 1 ) n 

~ ~ 2 1'' 

where  t=2n,w  = w(i),v  = w(i),  and  ' = 

n n n n ’ n 1 + ^/n 


Kquation  (B.16)  Koverninu  plastie  motion  has  tlio  solutions  lor 


I)  t I < 2 n; 


V + w .J./2 
, - , ' n n ' 

w = w cosv  (t  - t ) 4-  sins^  (1  - t)  e 

n n n V n 

n 


2 


P 1 o 
- p^^(-)e 

n J 


,(t-r)/2 


sin^J  ( t - t)(1t 
n 


(B. 19) 


, ( n - 1 ) ( k k,  n - w ) - . n/2  > (>.  2 

n 2 IB  o 


V . w 4 2) 

- , ' n n ' 

w = w eoshv  (t  - t)  . — sinliv  (t  - t)  e 

n n n n 

n 


-Li(t-t)  2 


- ( t — ) 2 

■)<'  sinhv  (I  - ')d- 


r 2 2 

wlierc  V = (p,/2)  - ' 


2 ■>  - 2 '■ 
X = k (n“  - l)(k  k n - w ) - an,  2 ^ 0 , X - (■,,  2) 

n 2 1 3 o n 


2 2 

or  V = X + (u./2) 

n n 


2 r 
\ = - 
n 


2 2 - 

k (n  - l)(k  k,n  - w ) - ;,n/2  > 0 

2 1 3 o 


Equation  (B.16)  f^overnin^  plastic  motion  has  the  solutions  lor 
0 < t < 2/n  < t or  0 < 2 n < ( •'  t: 


* * n 

w cos>  (t  “ t ) ^ - — sin\  (t 
n n \ n 

n 


(1  - \x/ri) 


r \ ■>:: 


( T ) s i n ( t 
n 


( B . 2 1) 


2 2 2 - 

whore  ' = k(n  - l)(kkn  - w)/(l.  n)  n 

n 2 13  o 


w = w cos''  (t  - t ) ► “ — sinhX  (t  - I ) 

n n n 1 n 


^ 1."" 


( - ) s i nh''  (1  - ')  ti- 
ll 


(B.22) 
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